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ABSTRACT 

We present a systematic analysis of the spectral and temporal properties of 17 
gamma-ray bursts (GRBs) co-detected by Gamma-Ray Monitor (GBM) and Large Area 
Telescope (LAT) on board the Fermi satellite by May 2010. We performed a time- 
resolved spectral analysis of all the bursts with the finest temporal resolution allowed 
by statistics, in order to reduce temporal smearing of different spectral components. 
We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with 
the classical "Band" function over the entire Fermi spectral range, which may sug- 
gest a common origin for emissions detected by LAT and GBM. GRB 090902B and 
GRB 090510 require the superposition of a MeV component and an extra power law 
component, with the former having a sharp cutoff above E p . For GRB 090902B, this 
MeV component becomes progressively narrower as the time bin gets smaller, and can 
be fit with a Planck function as the time bin becomes small enough. In general, we 
speculate that phenomenologically there may be three elemental spectral components 
that shape the time-resolved GRB spectra: a Band-function component (e.g. in GRB 
080916C) that extends in a wide energy range and does not narrow with decreasing time 
bins, which may be of non-thermal origin; a quasi-thermal component (e.g. in GRB 
090902B) with the spectra progressively narrowing with reducing time bins; and another 
non-thermal power law component extending to high energies. The spectra of differ- 
ent bursts may be decomposed into one or more of these elemental components. We 
compare this sample with the BATSE sample and investigate some correlations among 
spectral parameters. We discuss the physical implications of the data analysis results 
for GRB prompt emission, including jet composition (matter-dominated vs. Poynting- 
flux-dominated outflow), emission sites (internal shock, external shock or photosphere), 
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as well as radiation mechanisms (synchrotron, synchrotron self-Compton, or thermal 
Compton upscattering) . 

Subject headings: gamma-ray burst: general 

1. Introduction 

Although observationally accessed much earlier, GRB prompt emission is stil less understood 
than afterglow. The fundamental uncertainties lie in the following three poorly known important 
properties of GRBs (e.g. Zhang & Meszaros 2004 for a review): (1) Ejecta composition: Are the 
ejecta mostly composed of baryonic matter or a Poynting flux? (2) Energy dissipation site: Is the 
emission from internal shocks (Rees & Meszaros 1994; Kobayashi et al. 1997), the photosphere 
(Paczynski 1986; Goodman 1986; Meszaros & Rees 2000; Meszaros et al. 2002; Pe'er 2008), some 
magnetic dissipation regions (Lyutikov & Blandford 2003; Zhang & Yan 2011), or the external 
shock (Rees k, Meszaros 1992; Meszaros & Rees 1993; Dermer & Mitman 1999)? (3) Is the radiation 
mechanism synchrotron/jitter radiation (Meszaros et al. 1994; Medvedev 2000), synchrotron self- 
Compton (Panaitescu & Meszaros 2000; Kumar & McMahon 2008), or Comptonization of thermal 
photons (e.g. Thompson 1994; Pe'er et al. 2005, 2006; Beloborodov 2010; Lazzati & Begelman 
2009)? 

Before Fermi, understanding of GRB prompt emission has progressed slowly. Observations 
of early X-ray afterglows by Swift revealed a steep decay phase that is smoothly connected to 
prompt emission (Tagliaferri et al. 2005; Barthelmy et al. 2005), which suggests that the prompt 
emission region is detached from the afterglow emission region, and that the prompt emission 
site is "internal" (Zhang et al. 2006). Other than this, the properties of prompt emission were 
poorly constrained. The main factor that hampers progress has been the narrow energy band of 
the gamma-ray detectors of previous missions. Theoretical models usually predict rich features 
in the prompt spectra (e.g. Pe'er et al. 2006, see Zhang 2007; Fan & Piran 2008 for reviews on 
high energy emission processes). However, within the narrow observational spectral window, these 
features cannot be fully displayed. Instead, most previous spectral analyses revealed an empirical 
"Band" -function (Band et al. 1993), which is a smoothly-joint broken power law, whose physical 
origin is not identified. For the bright BATSE GRB sample, the typical low and high energy photon 
indices are distributed around a ~ —1 and (3 ~ —2.2, respectively, while the spectral peak energy 
E p is distributed around 200-300 keV (Preece et al. 2000). Later observations suggested that the 
distribution of E p can be much wider, extending to a few keV in the soft regime for X-ray flashes 
(Sakamoto et al. 2005) and to greater than 100 MeV in the hard regime (e.g. > 170 MeV for 
GRB 930506, Kaneko et al. 2008). Some BATSE GRBs were also detected by EGRET in the GeV 
range (Kaneko et al. 2008). For example, It was found that the GeV emission can last much longer 
than the MeV emission (e.g. GRB 940217, Hurley et al. 1994), and that it can form a distinct 
spectral component (e.g. GRB 941017, Gonzalez et al. 2003). In the softer regime, an X-ray excess 
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component with respect to the Band function was discovered in some BATSE GRBs (Preece et al. 
1996). However, the previous data were not adequate to place meaningful constraints on the three 
main questions discussed above. 

The Fermi satellite ushered in a new era of studying GRB prompt emission. The two instru- 
ments on board Fermi, the Gamma-ray Burst Monitor (GBM; Meegan et al. 2009) and the Large 
Area Telescope (LAT; Atwood et al. 2009), provide an unprecedented spectral coverage for 7 orders 
of magnitude in energy (from ~8 keV to ~300 GeV). Since the beginning of GBM/LAT science 
operation in August 2008 and as of the writing of this paper (May 2010), there have been 17 GRBs 
co-detected by LAT and GBM, with a detection rate comparable to the expectation assuming that 
the LAT-band emission is the simple extrapolation of the Band spectrum to the GeV range (Abdo 
et al. 2008; Lii et al. 2010). As will be shown below, the Band- function fits apply to most LAT 
GRBs, although some outliers do exist. Broad band spectral analyses have been published by the 
Fermi team for several individual GRBs, e.g. GRB 080916C (Abdo et al. 2009a), GRB 090510 
(Abdo et al. 2009b, Ackermann et al. 2010), GRB 090902B (Abdo et al. 2009c, Ryde et al. 2010), 
GRB 080825C (Abdo et al. 2009d), and GRB 081024B (Abdo et al. 2010a), which revealed several 
interesting features, such as the nearly featureless Band spectra covering 6 orders of magnitude in 
all epochs for GRB 080916C, the existence of an extra power law component extending to high 
energies in GRB 090510 and GRB 090902B, the existence of a quasi-thermal emission component 
in GRB 090902B, the delayed onset of the LAT-band emission with respect to the GBM-band 
emission, as well as an extended rapidly decaying GeV afterglow for most GRBs. 

These discoveries have triggered a burst of theoretical investigations of GRB prompt emission. 
Zhang Sz Pe'er (2009) argued that the lack of a thermal component in the nearly featureless spectra 
of GRB 080916C suggests a Poynting flux dominated flow for this burst. The conclusion was 
strengthened by a follow up study of Fan (2010, see also Gao et al. 2009). On the other hand, 
the quasi-thermal component in GRB 090902B (Ryde et al. 2010) is well-consistent with the 
photosphere emission of a hot fireball (Pe'er et al. 2010, Mizuta et al. 2010), suggesting that 
the burst is not highly magnetized. The possibility that the entire Band function spectrum is 
photosphere emission was discussed by several authors (Fan 2009; Toma et al. 2010; Beloborodov 
2010; Lazzati & Begelman 2010; Ioka 2010). These models have specific predictions that can be 
tested by the available data. In the high energy regime, Kumar & Barniol Duran (2009, 2010), 
Ghisellini et al. (2010) and Wang et al. (2010) suggested that the GeV afterglow is of external shock 
origin, which requires some unconventional parameters (Li 2010a; Piran & Nakar 2010). On the 
other hand, the fact that LAT emission is the natural spectral extension of GBM emission in some 
GRBs suggests that the GeV emission may be of an internal origin similar to MeV emission (Zhang 
& Pe'er 2009). Finally, the delayed onset of the GeV emission has been interpreted as emergence of 
the upscattered cocoon emission (Toma et al. 2009), synchrotron emission from shock accelerated 
protons (Razzaque et al. 2009), as well as delayed residual internal shock emission (Li 2010b). 
Again these models have specific predictions that may be tested by a detailed analysis of the data. 

Our goal is to systematically analyze the GRB data collected by the Fermi mission, aiming 
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at addressing some of the above mentioned problems in prompt GRB emission physics. Here we 
report the first paper in the series, which focuses on a comprehensive analysis of the GRBs that 
were co-detected by LAT and GBM. This sample has a much broader spectral coverage than the 
GBM-only GRBs, and therefore carries much more information about GRB prompt emission. The 
plan of the paper is the following. In Section 2 we describe the details of our sample selection 
and data analysis method. The data analysis results are presented in Section 3, with emphases 
on the unique features of some GRBs. We also present spectral parameter distributions and some 
possible correlations. In Section 4, we summarize the results and speculate on the existence of at 
least three elemental spectral components, and discuss their possible physical origins and possible 
combinations. In Section 5, we present the comparison between the emissions detected in the GBM 
band and that detected in the LAT band and discuss their physical connections. Our conclusions 
are summarized in Section 6 with some discussion. 



2. Sample and Data Reduction 

As of May 2010, 17 GRBs have been co-detected by Fermi LAT and GBM. Our sample includes 
all 17 GRBs (Table 1). We downloaded the GBM and LAT data for these GRBs from the public 
science support center at the official Fermi web site http://fermi.gsfc.nasa.gov/ssc/data/. An IDL 
code was developed to extract the energy-dependent lightcurves and time-dependent spectra for 
each GRB. This code was based on the Fermi RMFIT package (V3.3), the Fermi Science Tools 
(v9rl5p2) and the HEASOFT tools, which allows a computer to extract lightcurves and spectra 
automatically. The human involvement is introduced later to refine the analysis when needed. The 
code automatically performs the following tasks. 

1. Extract the background spectrum and lightcurve of the GBM data. Fermi records GBM 
data in several formats. For background reduction we use the CSPEC format data because 
it has a wider temporal coverage than the event data (time-tagged event, TTE, format). 
The background spectrum and lightcurve are extracted from some appropriate time intervals 
before and after the burs10, and the energy-dependent background lightcurves are modeled 
with a polynomial function B(E c h,t), where E c ^ is a specified energy band. 

2. Extract the source spectrum and lightcurve of the GBM data. This is done with the event 
(TTE) data. GBM has 12 Nal detectors (8 keV-1 MeV) and 2 BGO detectors (200 keV-40 
MeV). The overall signal-to-noise ratio (SNR) and peak count rate are calculated for each 
detector. The brightest Nal and BGO detectors are usually used for the analyses. If several 



An appropriate background time interval is typically when the lightcurve is "flat" with Poisson noise photons. 
For each burst, we select background time intervals as [-tb.i ,-ib,i] before the burst and [^,3,^,4] after the burst, where 
tb's are typically in the order of tens to hundreds of seconds. The exact values vary for different bursts due to their 
different brightnesses and the corresponding orbit slewing phases. 
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detectors have comparable brightnesses, all of them (usually 2-4 detectors) are taken for the 
analyses. By subtracting the background spectrum and lightcurve obtained in the previous 
step, the time-dependent spectra and energy-dependent lightcurves of the source in the GBM 
band are then obtained. 

3. Estimate the LAT-band background. Since only a small number of photons are detected by 
LAT for most GRBs, the background estimation should be performed cautiously. It is not 
straightforward to estimate an accurate LAT background using off-source regions around the 
trigger time. In our analyses, the LAT background is extracted using on-source region data 
long after the GBM trigger when the photon counts merge into a Poisson noise. 

4. Extract the LAT-band spectrum and lightcurve. Both "diffuse" and "transient" photons 
(level 0-3) are included. Since the LAT point spread function (PSF) strongly depends on the 
incident energy and the convention point of the tracker (Ohno et al. 2010), the photons are 
grouped into FRONT and BACK classes and their spectra are extracted separately based on 
different detector response files. The region of interest (ROI) that contains significant counts 
of LAT photons is further refined when necessary (Atwood et al. 2009; Abdo et al. 2009d). 

5. Extract the background-subtracted GBM and LAT lightcurves for different energy bands. In 
our analysis, the lightcurves are extracted in the following energy bands: 8-150 keV, 150-300 
keV, 300 keV-MeV, 1-30 MeV, and the LAT band (above 100 MeV). 

6. Make dynamically time-dependent spectral fits. Initially, the burst duration is divided in 
an arbitrary number of slices. The code then automatically refines the number of slices and 
the time interval for each slice, so that the photon counts in each bin (typically minimum 
20 counts for GBM spectra) give adequate statistics for spectral fitting (the reduced x 2 
is typically in the range of 0.75 - 1.5, a special case is GRB 090510, see Sect l3.2p . The 
time slices are defined to be be as small as possible as long as the extracted spectra satisfy 
these statistical criteria. The GBM spectra of the selected Nal and BGO detectors and the 
LAT "FRONT" and "BACK" type spectra are all extracted for each slice. These spectra, 
together with the corresponding response files (using the same one as the CSPEC data for 
LAT, or generated using gtrsp for GBM) are input into XSPEC (V 12.5.1) simultaneously to 
perform spectral fitting. The following spectral functions are considered (in order of increasing 
free parameters): single power law (PL), blackbody (BB, Planck function), power-law with 
exponential cutoff (CPL), and Band function. The models are tested based on the following 
principles: (1) If a one-component model can adequately describe the data (giving reasonable 
reduced x 2 , say, between 0.75 and 1.5), two-component models are not considered; (2) for one- 
component models, if a function with less free parameters can describe the data adequately, 
it is favored over the models with more parameters. (3) In addition, the Akaike's Information 
Criterion^! (AIC, Akaike 1974) is calculated to evaluate each model by considering both the 



2 AIC is defined by AIC = n In ( — ) + 2k, where n is the number of data points, k is the number of free parameters 
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fitting goodness (x 2 ) and the complexity of the model. We confirmed that the model with 
minimal AIC is the preferred model we choose based on the first two criteria. Nonetheless, 
since most GRBs have a Band- function spectra (see below), we also apply the Band function 
to those time bins that do not demand it in order to compare the fitting results between the 
Band function and other functions with less parameters (e.g. power law, blackbody, or power 
law with exponential cutoff). 

To assess the quality of a spectral fit, we use the traditional x 2 statistics. Due to the low count 
rate of LAT photons, we use the Gehrels (1986) weighting method in the high energy regime. We 
also employed the C-stat method (as used by the Fermi team), and found that the two methods 
usually give consistent results. We chose the x 2 method since it gives more reliable error estimates. 
All the model fitting parameters and x 2 statistics are presented in Table 2. For each burst, we 
present the time-dependent spectral parameters in the designated time bins defined by the statistics 
of spectral fitting, as well as the time-integrated spectral fit during the entire burst in the last row. 

3. Data Analysis Results 

The data analysis results are presented in Figs. [DI171 Each figure corresponds to one burst, 
and contains 10-11 panels. In the left panels, the lightcurves in 5 energy bands (8-150 keV, 150- 
300 keV, 300 keV-1 MeV, 1-30 MeV, and > 100 MeV) are presented in linear scale, together with 
the temporal evolution of the spectral parameters (a, j3, E p for Band function, kT for blackbody 
function, and T for single power law photon index). The top right panel is an example photon 
spectrum with model fitting, typically taken at the brightest time bin. The time-dependent model 
spectra are presented in the mid-right panel. The time-slices for the time-resolved spectral fitting 
are marked with vertical lines in the left panel lightcurves. In the bottom right panel, the GBM 
and LAT lightcurves are presented and compared in logarithmic scale. 

In the following, we discuss the results of several individual bright GRBs (Sect J3. Ill3.4p . and 
then discuss other GRBs in general (Sect l3.5|) . We then present statistics of spectral parameters 
(Sect J3.6|) and some possible correlations (Sect l3.7|) . 

3.1. GRB 080916C 

As shown in Fig|2j GRB 080916C is a long GRB with a duration ~ 66 s. The entire lightcurve 
can be divided into 6 segments. The smallest time bins during the brightest epochs (first two) are 
3.7 s and 5.4 s, respectively. This corresponds to a rest-frame time interval < 1 s (given its redshift 
4.35, Greiner et al. 2009). In all the time intervals, we found that the Band- function gives excellent 



of a particular model, and \ 2 ls the residual sum of squares from the estimated model (e.g. Shirasaki et al. 2008). 
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fits to the data, consistent with Abdo et al. (2009a). Initially there is a spectral evolution where 
the spectra "widen" with time (a hardening and (3 softening), but later the spectral parameters 
essentially do not evolve any more. We note that the steep (3 in the first time bin is mostly because 
of the non-detection in the LAT band. The tight upper limit above 100 MeV constrains the range 
of f3 not to be too hard. On the other hand, with GBM data alone, the data of in first time bin 

can be still fit as a Band function, with /3 2. similar to the values at later epochs. This 

suggests an alternative interpretation to the data: The high energy spectral index may be similar 
throughout the burst. The delayed onset of LAT-band emission may be because initially there is a 
spectral cutoff around 100 MeV, which later moves to much higher energies (e.g. above 13.2 GeV 
in the second time bin). 

It is interesting to note that the time integrated spectrum of GRB 080916C throughout the 
burst is also well fit with a Band function, where the spectral indices do not vary with time 
resolution. As an example, we present in Fig|18l the vF v spectra of GRB 080916C for three time 
bins with varying time resolution. Remarkably, the parameters do not vary significantly: a ~ —1.12, 

(3 ~ -2.25 for 3.5-8 s; a 1.0, (3 2.29 for 2-10 s; a 1.0, (3 2.27 for 0-20 s. This is in 

stark contrast with GRB 090902B discussed below. 



3.2. GRB 090510 

The short GRB 090510 was triggered with a precursor 0.5 s prior to the main burst. Two LAT 
photons were detected before the main burst. During the first time slice (0.45-0.5 s), no LAT band 
emission is detected, and the GBM spectrum can be well fit with a PL with an exponential cutoff 
(CPL hereafter). In the subsequent time slices, an additional PL component shows up, and the 
time-resolved spectra are best fit by the CPL + PL model. If one uses a Band + PL model to fit 
the data, the high energy spectral index (3 of the Band component cannot be constrained. If one 
fixes (3 to a particular value, it must be steeper than -3.5 in order to be consistent with the data. 
The CPL invoked in these fits has a low energy photon index Tcpl ~ —(0.6 — 0.8), which is very 
different from the case of a BB (where Tcpl ~ +1). On the other hand, the high-energy regime 
(exponential cutoff) is very similar to the behavior of a BB. 

Since this is a short GRB, we do not have enough photons to perform very detailed time- 
resolved spectral analysis. However, in order to investigate spectral evolution and the interplay 
between the MeV component and the extra PL component, we nonetheless make 4 time bins (see 
also Ackermann et al. 2010). As a result, the reduced x 2 °f each segment is outside the range of 
0.75 < x 2 /dof < 1.5 as is required for other GRBs. Our reduction results are generally consistent 
with those of the Fermi team (Abdo et al. 2009b; Ackermann et al. 2010). 
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3.3. GRB 090902B 

The spectrum of GRB 090902B is peculiar. Abdo et al. (2009c) reported that both the time- 
integrated and time-resolved spectra of this GRB can be fit with the Band+PL model. Ryde et al. 
(2010) found that the time-resolved spectra can be fit with a PL plus a multi-color blackbody model. 
This raises the interesting possibility that a blackbody-like emission component is a fundamental 
emission unit shaping the observed GRB spectra. 

In order to test this possibility, we carried out a series of time-resolved spectral analysis on the 
data (Fig flBj) , We first fit the time-integrated data within the time interval 0-20 s, and found that it 
can be fit with a model invoking a Band function and a power law, but with a poor x 2 /dof ~ 3.52. 
Compared with the Band component of other GRBs, this Band component is very narrow, with 
a ~ —0.58, f3 ~ —3.32. A CPL + PL model can give comparable fit, with Tcpl ~ —0.59. Next we 
zoom into the time interval 8.5 - 11.5 s, and perform spectral fits. The Band+PL and CPL+PL 
models can now both give acceptable fits, with parameters suggesting a narrower spectrum. For 
the Band+PL model, one has a ~ -0.07, /? ~ -3.69 with x 2 /dof = L26 - For tne CPL+PL model, 
one has Tcpl ~ —0.08 with x 2 /dof = 1.30. Finally we room in the smallest time bin (9.5 - 10 s) 
in which the photon counts are just enough to perform adequate spectral fits. We find that the 
Band + PL model can no longer constrain (3. The spectrum becomes even narrower, with a ~ 0.07 
and f3 < —5. The CPL+PL model can fit the data with a range of allowed Tcpl- In particular, 
if one fixes Tcpl ~ +1 (the Rayleigh- Jeans slope of a blackbody), one gets a reasonable fit with 
X 2 /dof = 0.92. This encourages us to suspect that a blackbody (BB) + PL model can also fit the 
data. We test it and indeed found that the model can fit the data with x 2 /dof = 1.11. These 
different models require different Tpl for the extra PL component, but given the low photon count 
rate at high energies, all these models are statistically allowed. Since the BB + PL model has less 
parameters than the CPL + PL and Band + PL models, we take this model as the simplest model 
for this smallest time interval. 

Next, we tried to divide the lightcurve of GRB 090902B into as many as time bins as possible 
so that the photon numbers in each time bin are large enough for statistically meaningful fits to 
be performed. Thanks to its high flux, we managed to divide the whole data set (0-30 s) into 62 
time bins. We find that the data in each time bin can be well fit by a BB+PL model, and that 
the BB temperature evolves with time. The fitting results are presented in Table 2 and FigllOl 
The time-integrated spectrum, however, cannot be fit with such a model (x 2 /dof = 14732/276). 
A Band+PL model gives a much improved fit, although the fit is still not statistically acceptable 
( x 2 /dof = 2024/275). The best fitting parameters are a = -0.83, /3 = -3.68, E p = 847 keV, and 
r = —1.85. Notice that the high energy photon index of the time-integrated Band spectrum is 
much steeper/softer than that observed in typical GRBs (Fig ]T9]) . 

In Fig l20l we display the lightcurves of both the thermal and the power-law components. It is 
found that the two components in general track each other. This suggests that the physical origins 
of the two components are related to each other (see Section 5 for more discussion). 
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An important inference from the analysis of GRB 090902B is that a Band-like spectrum can be 
a result of temporal superposition of many blackbody-like components. This raises the interesting 
possibility of whether all "Band" function spectra are superposed thermal spectra. From the 
comparison between GRB 090902B and GRB 080916C (FigJTSj). we find that such speculation is 
far-fetched. As discussed above, GRB 080916C shows no evidence of "narrowing" as the time bin 
becomes small (~ 1 s in the rest frame). In the case of GRB 090902B, a clear "narrowing" feature 
is seen. For the time integrated spectrum, GRB 080916C has a wide Band function (with a ~ —1.0, 
P ~ -2.27), while GRB 090902B (0-20 s) has a narrow Band function (with a = -0.58, /3 = -3.32) 
with worse reduced x 2 ■ Another difference between GRB 090902B and GRB 080916C is that the 
former has a PL component, which leverages the BB spectrum on both the low-energy and the 
high-energy ends to make a BB spectrum look more similar to a (narrow) Band function. GRB 
080916C does not have such a component, and the Band component covers the entire Fermi energy 
range (GBM &; LAT). We therefore conclude that GRB 090902B is a special case, whose spectrum 
may have a different origin from GRB 080916C (and probably most other LAT GRBs as well, see 
Sect 133] below). 



3.4. GRB 090926A 

This is another bright long GRB with a duration ~ 20 s. In our analysis, the lightcurve is 
divided into 9 segments. The Band function gives an acceptable fit to all the time bins (Fig llip . 
We however notice a flattening of f3 after ~ 11 s after the trigger. Also the Band function fit gives 
a worse reduced x 2 (although still acceptable) after this epoch. Since our data analysis strategy is 
to go for the simplest models, we do not explore more complicated models that invoke Band + PL 
or Band + CPL (as is done by the Fermi team, Abdo et al. 2010b). In any case, our analysis does 
not disfavor the possibility that a new spectral component emerges after ~ 11 s since the trigger 
(Abdo et al. 2010b). 



3.5. Other GRBs 

The time resolved spectra of other 13 GRBs are all adequately described by the Band function, 
similar to GRB 080916C. The Band-function spectral parameters are generally similar to GRB 
080916C. It is likely that these GRBs join GRB 080916C forming a "Band-only" type GRBs. In 
the current sample of 17 GRBs, only GRB 090510, GRB 090902B and probably GRB 090926A 
do not belong to this category and have an extra PL component extending to high energies. One 
caveat is that some GRBs in the sample are not very bright, so that we only managed to divide 
the lightcurves into a small number of time bins (e.g. 3 bins for GRB 080825C, 1 bin for GRB 
081024B, 3 bins for GRB 090328, 3 bins for GRB 091031, 2 bins for GRB 100225A, and 1 bin for 
GRB 100325A). So one cannot disfavor the possibility that the observed spectra are superposition 
of narrower components (similar to GRB 090902B). However, at comparable time resolution, GRB 
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090902B already shows features that are different from these GRBs: (1) the Band component is 
"narrower", and (2) there is an extra PL component. These two features are not present in other 
GRBs. We therefore suggest that most LAT/GBM GRBs are similar to GRB 080916C. 

3.6. Spectral parameter distributions 

Since the time-resolved spectra of most GRBs in our sample can be adequately described as a 
Band function, we present the distributions of the Band function parameters in this section. Since 
their MeV component may be of a different origin, GRB 090510 and GRB 090902B are not included 
in the analysis. 

The distributions of the spectral parameters a, f3, and E p are presented in Figure \19\ with a 
comparison with those of the bright BATSE GRB sample (Preece et al. 2000). It is found that 
the distributions peak at a = —0.9, /3 = —2.6, and E p ~ 781 keV, respectively. The a and j3 
distributions are roughly consistent with those found in the bright BATSE GRB sample (Preece 
et al. 2000). The E p distribution of the current sample has a slightly higher peak than the bright 
BATSE sample (Preece et al. 2000). This is likely due to a selection effect, namely, a higher E p 
would favor GeV detections. 

3.7. Spectral parameter correlations 

For time-integrated spectra, it was found that E p is positively correlated with the isotropic 
gamma-ray energy and the isotropic peak gamma-ray luminosity (Amati et al. 2002; Wei & Gao 
2003; Yonetoku et al. 2004). For time resolved spectra, E p was also found to be generally correlated 
with flux (and therefore luminosity, Liang et al. 2004), although in individual pulses, both a 
decreasing E p pattern and a .Ep-tracking-flux pattern have been identified (Ford et al. 1995; Liang 
& Kargatis 1996; Kaneko et al. 2006; Lu et al. 2010). 

In Fig|2Tl we present the E^-luminosity relations. Fig|21~a is for the global E p — L p iso correla- 
tion. Seven GRBs in our sample that have redshift information (and hence, the peak luminosity) 
are plotted against previous GRBs (a sample presented in Zhang et al. 2009). Since the correlation 
has a large scatter, all the GBM/LAT GRBs follow the same correlation trend. In particular, GRB 
090902B, whose E p is defined by the BB component, also follows a similar trend. This suggests that 
even if there may be two different physical mechanisms to define a GRB's E p , both mechanisms 
seem to lead to a broad E p — L p iso relation. It is interesting to note that the short GRB 090510 
(the top yellow point), even located at the upper boundary of the correlation, is still not an outlier. 
This is consistent with the finding (Zhang et al. 2009; Ghirlanda et al. 2009) that long/short GRBs 
are not clearly distinguished in the L p iso — E p domain. 

In Fig]21b. we present the internal E p — Lry^ so correlation. It is interesting to note that although 
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with scatter, the general positive correlation between E p and L 7) i SO as discovered by Liang et al. 
(2004) clearly stands. More interestingly, the BB-defined E p (in GRB 090902B) follows a similar 
trend to the Band-defined E p (e.g. in GRB 080916C and GRB 090926A), although different bursts 
occupy a different space region in the E p — L 7i ; so plane. 

In Fig l22l we present various pairs of spectral parameters in an effort to search for possible 
new correlations. The GRBs with redshift measurements are marked in colors, while those without 
redshifts are marked in gray with an assumed redshift z = 1. In order to show the trend of evolution, 
points for same burst are connected, with the beginning of evolution marked as a circle. 

No clear correlation pattern is seen in the E p — a and E p — (3 plots. Interestingly, a preliminary 
trend of correlation is found in the following two domains. 

• An a — f3 anti-correlation: Figl22a shows a rough anti-correlation between a and f3 in indi- 
vidual GRBs. This suggests that a harder a corresponds to a softer /3, suggesting a narrower 
Band function. In the time domain, there is evidence in some GRBs (e.g. GRB 080916C, 
GRB 090926A, and GRB 100414A, see Figs l2ll 11171) that the Band function "opens up" as 
time goes by, but the opposite trend is also seen in some GRBs (e.g. GRB 091031, Fig fT3j) . 
The linear Pearson correlation coefficients for individual bursts are insert in Fig]22a 

• A flux-a correlation: Fig|22b shows a rough correlation between flux and a. Within the 
same burst, there is rough trend that as the flux increases, a becomes harder. The linear 
Pearson correlation coefficients for individual bursts are presented in Figl22b inset. One 
possible observational bias is that when flux is higher, one tends to get a smaller time slice 
based on the minimum spectral analysis criterion. If the time smearing effect can broaden the 
spectrum, then a smaller time slice tends to give a narrower spectrum, and hence, a harder 
a. This would be relevant to bursts similar to GRB 090902B, but not bursts similar to GRB 
080916C (which does not show spectral evolution as the time resolution becomes finer). More 
detailed analyses of bright GRBs can confirm whether such a correlation is intrinsic or due 
to the time resolution effect discussed above. 

Several caveats should be noted for these preliminary correlations: First, some bursts do not 
obey these correlations, so the correlations, if any, are not universal; Second, the currently chosen 
time bins are based on the requirement for adequate spectral analyses, so the time resolution varies 
in different bursts. For some bright bursts, a burst pulse can be divided into several time bins, 
while in some faint others, a time bin corresponds to the entire pulse; Third, the current sample 
is still too small. A time-resolved spectral analysis for more bright GBM GRBs may confirm or 
dispute these correlations. 
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4. Elemental Spectral Components and their physical origins 

4.1. Three phenomenologically identified elemental spectral components 

The goal of our time-resolved spectral analysis is to look for "elemental" emission units that 
shape the observed GRB prompt gamma-ray emission. In the past it has been known that time- 
integrated GRB spectra are mostly fit by the Band function (Band et al. 1993). However, whether 
this function is an elemental unit in the time-resolved spectra is not known. One speculation is 
that this function is the superposition of many simpler emission units. If such a superposition relies 
on adding the emission from many time slices, then these more elemental units should show up as 
the time bins become small enough. 

One interesting finding of our time-resolved spectral analyses is that the "Band" -like spectral 
component seen in GRB 090902B is different from that seen in GRB 080916C and some other 
Band-only GRBs. While the Band spectral indices of GRB 080916C essentially do not change 
as the time bins become progressively smaller, that of GRB 090902B indeed show the trend of 
"narrowing" as the time bin becomes progressively smaller. With the finest spectral resolution, 
GRB 090902B spectra can be fit by the superposition of a PL component and a CPL function, 
including a Planck function. Even for the time- integrated spectrum, the "Band" -like component in 
GRB 090902B appears "narrower" than that of GRB 080916C. All these suggest that the "Band"- 
like component of GRB 090902B is fundamentally different from that detected in GRB 080916C 
and probably also other Band-only GRBg^. Similarly, the time-resolved spectra of the short GRB 
090510 can be well fit by the superposition of a PL component and a CPL spectrum (although not 
a Planck function). The PL component extends to high energies with a positive slope in vF v . The 
CPL component may be modeled as a multi-color blackbody spectrum. We therefore speculate 
that the MeV component of GRB 090510 is analogous to that of GRB 090902B. 

Phenomenologically, the power law component detected in GRB 090902B and GRB 090510 
is an extra component besides the Band-like component. Such a component may have been also 
detected in the BATSE-EGRET burst GRB 941017 (Gonzalez et al. 2003), and may also exist in 
GRB 090926A at later epochs. 

We therefore speculate that phenomenologically there might be three elemental spectral compo- 
nents that shape the prompt gamma-ray spectrum. These include: (I) a Band function component 
("Band" in abbreviation) that covers a wide energy range (e.g. 6-7 orders of magnitude in GRB 
080916C) and persists as the time bins become progressively smaller. It shows up in GRB 080916C 
and 13 other LAT GRBs; (II) a quasi-thermal component ("BB" in abbreviatior0) which becomes 



3 Our finest time interval is around Is in the rest frame of the burst. Theoretically, how time-integrated spectra 
broaden with increasing time bins is subject further study. Our statement is therefore relevant for time resolution 
longer than Is. 

4 Notice that the abbreviation "BB" here not only denotes blackbody, but also includes various modifications to 
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progressively narrower as the time bin becomes smaller, and eventually can be represented as a 
blackbody (or multi-color blackbody) component as seen in GRB 090902B; (III) a power law com- 
ponent ("PL" in abbreviation) that extends to high energy as seen in GRBs 090902B and 090510, 
which has a positive slope in the vF v spectrum and should have an extra peak energy (E p ) at an 
even higher energy that is not well constrained by the data. 

Figure [23] is a cartoon picture of the vF v spectrum that includes all three phenomenologically 
identified elemental spectral components. The time resolved spectra of the current sample can be 
understood as being composed of one or more of these components. For example, GRB 080916C 
and other 13 GRBs have Component I (Band), GRB 090902B and probably GRB 090510 have 
Components II (BB) and III (PL), and GRB 0900926A has Component I initially, and may have 
components I and III at later times. 

4.2. Possible physical origins of the three spectral components 

4-2.1. Band Component 

The fact that the this component extends through a wide energy range (e.g. 6-7 orders of 
magnitude for GRB 080916C) strongly suggests that a certain non-thermal emission mechanism 
is in operation. This demands the existence of a population of power-law-distributed relativistic 
electrons, possibly accelerated in internal shocks or in regions with significant electron heating, 
e.g. magnetic dissipation. In the past there have been three model candidates for prompt GRB 
emission: synchrotron emission, synchrotron self-Compton (SSC), and Compton upscattering of a 
thermal photon source. In all these models the high energy PL component corresponds to emission 
from a PL-distributed electron population. The spectral peak energy E p may be related to the 
minimum energy of the injected electron population, an electron energy distribution break, or the 
peak of the thermal target photons. 

Most prompt emission modeling (Meszaros et al. 1994; Pilla & Loeb 1998; Pe'er & Waxman 
2004a; Razzaque et al. 2004; Pe'er et al. 2006; Gupta & Zhang 2007) suggest that the overall 
spectrum is curved, including multiple spectral components. Usually a synchrotron component 
is accompanied by a synchrotron self-Compton (SSC) component. For matter-dominated fireball 
models, one would expect the superposition of emissions from the photosphere and from the internal 
shock dissipation regions. As a result, the fact that 14/17 (~ 80%) of GRBs in our sample have a 
Band-only spectrum is intriguing. The three theoretically expected spectral features, i.e. the quasi- 
thermal photosphere emission, the SSC component (if the MeV component is of synchrotron origin), 
and a pair-production cutoff at high energies, are all not observed. This led to the suggestion that 
the outflows of these GRBs are Poynting flux dominated (Zhang & Pe'er 2009). Within such a 



the blackbody spectrum such as multi-color blackbody. 
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picture, the three missing features can be understood as the following: (1) Since most energy is 
carried in magnetic fields and not in photons, the photosphere emission (BB component) is greatly 
suppressed; (2) Since the magnetic energy density is higher than the photon energy density, the 
Compton Y parameter is smaller than unity, so that the SSC component is naturally suppressed; 
(3) A Poynting flux dominated model usually has a larger emission radius than the internal shock 
model (Lyutikov & Blandford 2003 for current instability and Zhang & Yan 2011 for collision- 
induced magnetic reconnection/turbulence). This reduces the two-photon annihilation opacity and 
increases the pair cutoff energy. This allows the Band component extend to very high energy (e.g. 
13.2 GeV for GRB 080916C). 

Another possibility, advocated by Beloborodov (2010) and Lazzati & Begelman (2010) in view 
of the Fermi data (see also discussion by Thompson 1994; Rees & Meszaros 2005; Pe'er et al. 2006; 
Giannios & Spruit 2007; Fan 2009; Toma et al. 2010; and Ioka 2010), is that the Band component is 
the emission from a dissipative photosphere. This model invokes relativistic electrons in the regions 
where Thomson optical depth is around unity, which upscatter photosphere thermal photons to 
high energies to produce a power law tail. This model can produce a Band-only spectrum, but has 
two specific limitations. First, the high energy power law component cannot extend to energies 
higher than GeV in the cosmological rest frame, since for effective upscattering, the emission region 
cannot be too far above the photosphere. The highest photon energy detected in GRB 080916C 
is 13.2 GeV (which has a rest-frame energy ~ 70 GeV for its redshift z = 4.35). This disfavors 
the dissipative photosphere model. This argument applies if the LAT-band photons are from the 
same emission region as the MeV photons, as suggested by the single Band function spectral fits. 
It has been suggested that the LAT emission during the prompt phase originates from a different 
emission region, e.g. the external shock (Kumar & Barniol Duran 2009; Ghisellini et al. 2010). 
This requires that the two distinct emission components conspire to form a nearly featureless Band 
spectrum in all temporal epochs, which is contrived. As will be shown in Sect J5.2l later, there is 
compelling evidence that the LAT emission during the prompt emission phase is of an internal 
origin. In particular, the peak of the GeV lightcurve of GRB 080916C coincides with the second 
(the brightest) peak of GBM emission, and the 13.2 GeV photon coincides with another GBM 
lightcurve peak. All these suggest an internal origin of the GeV emission during the prompt phase. 

The second limitation of the dissipative photosphere model is that the photon spectral index 
below E p is not easy to reproduce. The simplest blackbody model predicts a Rayleigh-Jeans spec- 
trum a = +1. By considering slow heating, this index can be modified as a = +0.4 (Beloborodov 
2010). Both are much harder from the observed a ~ — 1 value. In order to overcome this difficulty, 
one may appeal to the superposition effect, i.e. the observed Band spectrum is the superposition 
of many fundamental blackbody emission units (e.g. Blinnikov et al. 1999; Toma et al. 2010; 
Mizuta et al. 2010; Pe'er & Ryde 2010). However, no rigorous calculation has been performed to 
fully reproduce the a = — 1 spectrum. Pe'er &: Ryde (2010) show that when the central engine 
energy injection is over and the observed emission is dominated by the high-latitude emission, an 
a = — 1 can be reproduced with the flux decaying rapidly with oc t~ 2 . During the phase when the 
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central engine is still active, the observed emission is always dominated by the contribution along 
the line of sight, which should carry the hard low energy spectral index of the blackbody function. 
Observationally, the Band component spectral indices are not found to vary when the time bins 
are reduced (in stark contrast to the narrow Band-like component identified in GRB 090902B). 
This suggests that at least the temporal superposition of many blackbody radiation units is not 
the right interpretation for this component. 

4-2.2. Quasi- Thermal (BB) Component 

The MeV component in GRB 090902B narrows with reduced time resolution and eventually 
turns into being consistent with a blackbody (or multi-color blackbody) as the time bin becomes 
small enough. This suggests a thermal origin of this component. Within the GRB content, a 
natural source is the emission from the photosphere where the photons advected in the expanding 
relativistic outflow turn optically thin for Compton scattering. In fact, the original fireball model 
predicts a quasi-thermal spectrum (Paczynski 1986; Goodman 1986). In the fireball shock model, 
such a quasi-blackbody component is expected to be associated with the non-thermal emission 
components (Meszaros & Rees 2000; Meszaros et al. 2002; Daigne & Mochkovitch 2002; Pe'er et 
al. 2006). 

Some superposition effects may modify the thermal spectrum to be different from a pure Planck 
function. The first is the temporal smearing effect. If the time bin is large enough, one samples 
photosphere emission from many episodes, and hence, the observed spectrum should be a multi- 
color blackbody. This effect can be diminished by reducing the time bin for time-resolved spectral 
analyses. GRB 090902B is such an example. The second effect is inherited in emission physics of 
relativistic objects. At a certain epoch, the observer detects photons coming from different latitudes 
from the line of sight, with different Doppler boosting factors. The result is an intrinsic smearing 
of the Planck function spectrum. Pe'er & Ryde (2010) have shown that after the central engine 
activity ceases, the high-latitude emission effect would give an a ~ —1 at late times, with a rapidly 
decaying flux F u oc t~ 2 . This second superposition effect is intrinsic, and cannot be removed by 
reducing the time bins. 

The case of the thermal component is most evidenced in GRB 090902B, and probably also in 
GRB 090510. In both bursts, the MeV component can be well fit with a CPL + PL spectrum. 
The exponential cutoff at the high energy end is consistent with thermal emission with essentially 
no extra dissipation. For GRB 090902B, the low energy spectral index Tcpl is typically ~ 0, and 
can be adjusted to +1 (blackbody). For GRB 090510, Tcpl is softer (~ —0.7). Since it is a short 
GRB, the high-latitude effect may be more important. The softer low energy spectral index may 
be a result of the intrinsic high-latitude superposition effect (Pe'er & Ryde 2010). 
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4-2.3. Power-Law (PL) Component 

This component is detected in GRB 090902B and GRB 090510. Several noticeable properties 
of this component are: (1) For our small sample, this component is always accompanied by a low 
energy MeV component (likely the BB component). Its origin may be related to this low energy 
component; (2) It is demanded in both the low energy end and the high energy end, and amazingly 
the same spectral index can accommodate the demanded excesses in both ends. This suggests 
that either this PL component extends for 6-7 orders of magnitude in energy, or that multiple 
emission components that contribute to the excesses in both the low and high energy regimes have 
to coincide to mimic a single PL; (3) The spectral slope is positive in the vF v space, so that the 
main energy power output of this component is at even higher energies (possibly near or above 
the upper bound of the LAT band). 

Since the non-thermal GRB spectra are expected to be curved (Meszaros et al. 1994; Pilla & 
Loeb 1998; Pe'er & Waxman 2004a; Razzaque et al. 2004; Pe'er et al. 2006; Gupta & Zhang 2007; 
Asano & Terasawa 2009), the existence of the PL component is not straightforwardly expected. It 
demands coincidences of various spectral components to mimic a single PL component in the low and 
high energy ends. Pe'er et al. (2010) have presented a theoretical model of GRB 090902B. According 
to this model, the apparent PL observed in this burst is the combination of the synchrotron emission 
component (dominant at low energies), the SSC and Comptonization of the thermal photons (both 
dominant at high energies). A similar model was analytically discussed by Gao et al. (2009) within 
the context of GRB 090510. 

One interesting question is how Component III (PL) differs from Component I (Band). Since 
both components are non-thermal, they may not be fundamentally different. They can be two 
different manifestations of some non-thermal emission mechanisms (e.g. synchrotron and inverse 
Compton scattering) under different conditions. On the other hand, since Component III seems to 
be associated with Component II (BB) (e.g. in GRB 090902B and GRB 090510), its origin may 
be related to Component II. One possible scenario is that Component III (at least the part above 
component II) is the Compton-upscattered emission of Component II (e.g. Pe'er & Waxman 2004b 
for GRB 941017). The fact that the lightcurves of the BB component and the PL component of 
GRB 090902B roughly track each other (Fig l20p generally supports such a possibility. Within 
this interpretation, one must attribute the PL part below the thermal peak as due to a different 
origin (e.g. synchrotron, see Pe'er et al. 2010). Alternatively, Component I and III may be 
related to non-thermal emission from two different emission sites (e.g. internal vs. external or two 
different internal locations). Indeed, if the late spectra of GRB 090926A are the superposition of 
the components I and III, then both components can coexist, which may correspond to two different 
non-thermal emission processes and/or two different emission sites. 
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4.3. Possible Spectral Combinations of GRB Prompt Emission 

Using the combined GBM and LAT data, we have phenomenologically identified three ele- 
mental spectral components during the prompt GRB phase (Fig |23p . Physically they may have 
different origins (see above). One may speculate that all the GRB prompt emission spectra may 
be decomposed into one or more of these spectral components. It is therefore interesting to investi- 
gate how many combinations are in principle possible, how many have been discovered, how many 
should not exist and why, and how many should exist and remain to be discovered. We discuss the 
following possibilities in turn below (see Figl24l for illustrations). 

1. Component I (Band) only: 

This is the most common situation, which is observed in 14/17 GRBs in our sample exempli- 
fied by GRB 080916C. Either the BB and PL components do not exist, or they are too faint 
to be detected above the Band component. If the BB component is suppressed, these bursts 
may signify non-thermal emission from an Poynting flux dominated flow. 

2. Component II (BB) only: 

No such case exists in the current sample. GRB 090902B, and probably also GRB 090510, 
have a BB component, but it is accompanied by a PL component in both GRBs. It remains to 
be seen whether in the future a BB-only GRB will be detected, or whether a BB component is 
always accompanied by a PL component. Since the case of GRB 090902B is rare, we suspect 
that the BB-only GRBs are even rarer, if they exist at all. 

3. Component III (PL) only: 

Our PL component stands for the high energy spectral component seen in GRB 090902B 
and GRB 090510, which likely has a high E p near or above the boundary of the LAT band. 
Observationally, there is no solid evidence for such PL-only GRB^fl. In our current sample 
which covers the widest energy band, the PL component only exists in 2 out of 17 GRBs, 
and is found to be associated with the BB component. The luminosity of the PL component 
is found to roughly track that of the thermal component (Fig |20p . If the PL component is 
the Comptonization of a low energy photon source (e.g. the BB component), then PL-only 
GRBs may not exist in nature. 

4. I + II: 



5 Most of Swift GRBs can be fit with a PL (Sakamoto et al. 2008). However, this is due to the narrowness of the 
energy band of the gamma-ray detector BAT on board Swift. The E p of many Swift GRBs are expected to be located 
outside the instrument band. In fact, using a Band function model and considering the variation of E v within and 
outside the BAT band, one can reproduce the apparent hardness of Swift GRBs, and obtain an effective correlation 
between the BAT-band photon index and E p (e.g. Zhang et al. 2007; Sakamoto et al. 2009). If a GRB is observed 
in a wider energy band, the spectrum should be invariably curved. 
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Such a case is not found in our sample. If the Band component is the emission from the 
internal shocks and the BB component is the emission from the photosphere, then such a 
combination should exist and be common for fireball scenarios. An identification of such 
a case would confirm the non-thermal nature of the Band component (since the thermal 
component is manifested as the BB component). Observationally, an X-ray excess has been 
observed in 12 out of 86 (~ 14%) bright BATSE GRBs (Preece et al. 1996). This could 
be due to the contamination of a BB component in the X-ray regime. With the excellent 
spectral coverage of Fermi, we expect that such a spectral combination may be identified 
in some GRBs, even if technically it may be difficult because there are too many spectral 
parameters to constrain. 

5. I + III: 

Such a combination has not been firmly identified in our sample. Nonetheless, the spectral 
hardening of GRB 090926A after lis may be understood as the emergence of the PL com- 
ponent on top of the Band component seen before 11 s. Physically it may be related to two 
non-thermal spectral components or non-thermal emission from two different regions. 

6. II + III: 

Such a case is definitely identified in GRB 090902B, and likely in GRB 090510 as well. From 
the current sample, it seems that such a combination is not as common as the Band-only type, 
but nonetheless forms a new type of spectrum that deserves serious theoretical investigations. 
Physically, the high-energy PL component is likely the Compton up-scattered emission of the 
BB component, although other non-thermal processes (e.g. synchrotron and SSC) could also 
contribute to the observed emission (Pe'er et al. 2010). 

7. I + II + III: 

The full combination of all three spectral components (e.g. Fig ]23p is not seen from the current 
sample. In any case, in view of the above various combinations (including speculative ones), 
one may assume that the full combination of the three spectral components is in principle 
possible. Physically this may correspond to one photosphere emission component and two 
more non-thermal components (either two spectral components or non-thermal emission from 
two different regions). Nonetheless, technically there are too many parameters to constrain, 
so that identifying such a combination is difficult. 



5. LAT-band emission vs. GBM-band emission 

Besides the joint GBM/LAT spectral fits, one may also use temporal information to investigate 
the relationship between the emission detected in the GBM-band and that detected in the LAT 
band. In this section we discuss three topics: delayed onset of LAT emission, rough tracking 
behavior between GBM and LAT emissions, and long-lasting LAT afterglow. 
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5.1. Delayed onset of LAT emission 

The Fermi team has reported the delayed onset of LAT emission in several GRBs (GRBs 
080825C, 080916C, 090510, 090902B, Abdo et al. 2009a,b,c,d). Our analysis confirms all these 
results. In Table 1, we mark all the GRBs in our sample that show the onset delay feature. 

There have been several interpretations to the delayed onset of GeV emission discussed in the 
literature. Toma et al. (2009) suggested that GeV emission is the upscattered cocoon emission 
by the internal shock electrons. Razzaque et al. (2009) interpreted the GeV emission as the 
synchrotron emission of protons. Since it takes a longer time for protons to be accelerated and 
be cooled to emit GeV photons, the high energy emission is delayed. Li (2010b) interpreted GeV 
emission as the upscattered prompt emission photons by the residual internal shocks. 

Although it is difficult to test these models using the available data, our results give some 
observational constraints to these models. First, except GRBs 090510 and 090902B whose GeV 
emission is a distinct spectral component, other GRBs with onset delay still have a simple Band- 
function spectrum after the delayed onset. This suggests that for those models that invoke two 
different emission components to interpret the MeV and GeV components, one needs to interpret 
the coincidence that the GeV emission appears as the natural extension of the MeV emission to 
the high energy regime. 

For such delayed onsets whose GeV and MeV emissions form the same Band component, one 
may speculate two simpler explanations. One is that there might be a change in the particle 
acceleration conditions (e.g. magnetic configuration in the particle acceleration region). As shown 
in Sect J3.lt the early spectrum during the first time bin (before onset of LAT emission) of GRB 
080916C may be simply a consequence of changing the electron spectral index. One may speculate 
that early on the particle acceleration process may not be efficient, so that the electron energy 
spectral index is steep. After a while (the observed delay), the particle acceleration mechanism 
becomes more efficient, so that the particle spectral index reaches the regular value. The second 
possibility is that there might be a change in opacity. The GBM data alone during the first time bin 
gives a similar (3 as later epochs. It is possible that there might be a spectral cutoff slightly above 
the GBM band early on. A speculated physical picture would be that the particle acceleration 
conditions are similar throughout the burst duration, but early on the pair production opacity may 
be large (probably due to a lower Lorentz factor or a smaller emission radius), so that the LAT band 
emission is attenuated. The opacity later drops (probably due to the increase of Lorentz factor 
or the emission radius), so that the LAT band emission can escape from the GRB. Within such a 
scenario, one would expect to see a gradual increase of maximum photon energy as a function of 
time. Figure 25 shows the LAT photon arrival time distribution of GRB 080916C. Indeed one can 
see a rough trend of a gradual increase of the maximum energy with time. 

One last possibility is that the LAT band emission is dominated by the emission from the 
external shock, which is delayed with respect to the GBM-band prompt emission. This possibility 
is discussed in more detail below in Sec. 15.31 
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5.2. Rough tracking behavior 

Inspecting the multi-band lightcurves (FigsfTlfTTl left panels), for bright GRBs (e.g. 080916C, 
090217, 090323, 090902B) the LAT emission peaks seem to roughly track some peaks of the GBM 
emission (aside from the delayed onset for some of them). For example, the peak of the LAT 
lightcurve of GRB 080916C coincides with the second GBM peak. This is consistent with the 
spectral analysis showing that most time-resolved joint spectra are consistent with being the same 
(Band- function) spectral component. Even for GRB 090902B whose LAT band emission is from a 
different emission component from the MeV BB component, the emissions in the two bands also 
roughly track each other (Fig l20p . This suggests that the two physical mechanisms that power the 
two spectral components are related to each other. 

The rough tracking behavior is evidence against the proposal that the entire GeV emission 
is from the external forward shock (see Sec. 15.31 for more discussion). Within the forward shock 
model, the fluctuation in energy output from the central engine should be greatly smeared, since 
the observed flux change amplitude is related to AE/E <C 1 (where E is the total energy already 
in the balstwave, and AE is the newly injected energy from the central engine), rather than AE 
itself within the internal models. 



5.3. Long Term Emission in the LAT Band 

In order to study the long-term lightcurve behavior, we extract the GBM-band and LAT- 
band lightcurves in logarithmic scale and present them in the bottom right panel of FigsfTlfTTl We 
unevenly bin the LAT lightcurves with bin sizes defined by the requirement that the signal-to-noise 
ratio must be > 5. For a close comparison, we correspondingly re-bin the GBM lightcurves using 
the same bin sizes. Some GRBs (e.g. 080916C, 090510, 090902B, and 090926A) have enough 
photons to make a well sampled LAT lightcurve. 

In several GRBs, LAT emission lasts longer than GBM emission and decays as a single power 
law (Ghisellini et al. 2010). The decay indices of LAT emission are marked in the last panel of 
Figsm fTTl which can be also found in Table 3. Due to low photon numbers, it is impossible to carry 
out a time resolved spectral analysis. In any case, the LAT-band photon indices of long-term LAT 
emission are estimated and also presented in Table 3. In Table 1 we mark those GRBs with detected 
LAT emission longer than GBM emission and those without. The most prominent ones with long 
lasting LAT afterglow are GRBs 080916C, 090510, 090902B, and 090926A. Spectral analyses suggest 
that the LAT emission in GRBs 090510 and 090902B is a different spectral component from the 
MeV emission. The GBM lightcurves of these GRBs indeed follow a different trend by turning 
off sharply as compared with the extended PL decay in the LAT band. GRB 090926A, on the 
other hand, shows a similar decay trend in both GBM and LAT bands. GRB 080916C is special. 
Although the spectral analysis shows a single Band function component, the GBM lightcurve turns 
over sharply around 70-80 seconds, while the LAT emission keeps decaying with a single PL. 
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One caveat of LAT long-term lightcurves is that they depend on the level of background 
and time-bin selection. Due to the low count rate at late times, the background uncertainty can 
enormously change the flux level, and a different way of binning the data may change the shape 
of the lightcurve considerably. In our analysis, the background model is extracted from the time 
interval prior to the GBM trigger in the same sky region that contains the GRB. The bin-size is 
chosen to meet the 5a statistics to reduce the uncertainty caused by arbitrary binning. 

Our data analysis suggests a controversial picture regarding the origin of this GeV afterglow. 
Spectroscopically, the LAT-band emission is usually an extension of the GBM-band emission and 
forms a single Band-function component, suggesting a common physical origin with the GBM-band 
emission. If one focuses on the prompt emission lightcurves, the LAT-band activities seem to track 
the GBM-band activities. Even for GRB 090902B which shows a clear second spectral component, 
the PL component variability tracks that of the BB component well (Fig J2Q|) , suggesting a physical 
connection between the two spectral components. These facts tentatively suggest that at least 
during the prompt emission phase, the LAT-band emission is likely connected to the GBM-band 
emission, and may be of an "internal" origin similar to the GBM-band emission. 

It has been suggested that the entire GeV emission originates from the external shock (e.g. 
Kumar & Barniol Duran 2009a, 2009b; Ghisellini et al. 2010; Corsi et al. 2009). This idea is based 
on the power law temporal decay law that follows the prompt emission. Such a GeV afterglow 
scenario is not straightforwardly expected for the following reasons. First, before Fermi, afterglow 
modeling suggests that for typical afterglow parameters, the GeV afterglow is initially dominated 
by the synchrotron self-Compton component (Meszaros & Rees 1994; Dermer et al. 2000; Zhang & 
Meszaros 2001; Wei & Fan 2007; Gou & Meszaros 2007; Galli & Piro 2007; Yu et al. 2007; Fan et 
al. 2008), or by other IC processes invoking both forward and reverse shock electrons (Wang et al. 
2001). For very energetic GRBs such as GRB 080319B, one may expect a synchrotron-dominated 
afterglow all the way to an energy ~ 10 GeV (Zou et al. 2009; Fan et al. 2008). Second, the 
required parameters for the external shock are abnormal to interpret the data. For example, the 
magnetic field strength at the forward shock needs to be much smaller than equipartition, consistent 
with simply compressing the ISM magnetic field without shock amplification (Kumar & Barniol 
Duran 2010). This, in turn, causes a problem in accelerating electrons to a high enough energy 
to enable emission of GeV photons (Li 2010a; Piran & Nakar 2010). Moreover, the circumburst 
number density of these long GRBs are required to be much lower than that of a typical ISM 
(e.g., Kumar & Barniol Duran 2010), which challenges the collapsar model. Finally, observed GeV 
decay slope is typically steeper than the predictions invoking a standard adiabatic forward shock 
(e.g. Figs l2|8ll0lllll7"l see also Ghisellini et al. 2010). One needs to invoke a radiative blastwave 
(Ghisellini et al. 2010) or a Klein-Nishina cooling-dominated forward shock (Wang et al. 2010) to 
account for the steepness of the decay slope. 

The external shock model to interpret the entire GeV emission is challenged by the following 
two arguments. First, the GeV lightcurve peak coincides the second peak of the GBM lightcurve for 
GRB 080916C. This requires a fine-tuned bulk Lorentz factor of the fireball to make the deceleration 
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time coincide the epoch of the second central engine activity. This is highly contrived. Second, the 
external shock component should not have decayed steeply while the prompt emission is still on 
going. To examine this last point, we have applied the shell-blastwave code developed by Maxham 
& Zhang (2009) to model the blastwave evolution of GRB 080916C using the observed data by 
assuming that the outflow kinetic energy traces the observed gamma-ray lightcurve (assuming a 
constant radiation efficiency). The resulting LAT-band lightcurve always displays a shallow decay 
phase caused by refreshing the forward shock by materials ejected after the GeV lightcurve peak 
time even for a radiative blastwave, in stark contrast to the data. This casts doubts on the external 
shock origin of GeV emission during the prompt phase (Maxham et al. 2011). We note that detailed 
modeling of GRB 090510 (He et al. 2010) and GRB 090902B (Liu & Wang 2011) with the external 
shock model both suggests that the prompt GeV emission cannot be interpreted as the emission 
from the external forward shock. 

Collecting the observational evidence and the theoretical arguments presented above, we sug- 
gest that at least during the prompt emission phase (when GBM-band emission is still on), the 
LAT-band emission is not of external forward shock origin. 

After the GBM-band prompt emission is over, the LAT-band emission usually decays as a PL. 
We note that the long-term GeV lightcurve can be interpreted in more than one way. (1) If one 
accepts that the prompt GeV emission is of internal origin, one may argue that the external shock 
component sets in before the end of the prompt emission and thereafter dominates during the decay 
phase (Maxham et al. 2011). This requires arguing for coincidence of the same decaying index for 
the early internal and the late external shock emission. Considering a possible superposition effect 
(i.e. the observed flux during the transition epoch includes the contributions from both the internal 
and external shocks) , this model is no more contrived than the model that interprets prompt GeV 
emission as from external shocks, which requires coincidence of internal emission spectrum and 
the external shock emission spectrum to mimic the same Band spectrum in all time bins (Kumar 
& Barniol Duran 2009). (2) An alternative possibility is to appeal to an internal origin of the 
entire GeV long-lasting afterglow, which reflects the gradual "die-off" of the central engine activity. 
The difficulty of such a suggestion is that it must account for the different decaying behaviors 
between the GBM-band emission and LAT-band emission in some (but not all) GRBs (e.g. GRB 
080916C). To differentiate between these possibilities, one needs a bright GRB co-triggered by 
Fermi LAT/GBM and Swift BAT, so that an early Swift XRT lightcurve is available along with the 
early GeV lightcurve. The external-shock-origin GeV afterglow should be accompanied by a PL 
decaying early X-ray lightcurve (Liang et al. 2009) instead of the canonical steep-shallow-normal 
decaying pattern observed in most Swift GRBs (Zhang et al. 2006; Nousek et al. 2006; O'Brien et 
al. 2006). A violation of such a prediction would suggest an internal origin of the GeV afterglow. 
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6. Conclusions and discussion 

We have presented a comprehensive joint analysis of 17 GRBs co-detected by Fermi GBM 
and LAT. We carried out a time-resolved spectral analysis of all the bursts with the finest tem- 
poral resolution allowed by statistics, in order to reduce temporal smearing of different spectral 
components. Our data analysis results can be summarized as the following: 

• We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with the 
classical "Band" function over the entire Fermi spectral range, which may suggest a common 
origin for emissions detected by LAT and GBM. GRB 090902B and GRB 090510 are found 
to be special in that the data require the superposition between a MeV component and an 
extra power law component, and that the MeV component has a sharp cutoff above E p . More 
interestingly, the MeV component of GRB 090902B becomes progressively narrower as the 
time bin gets smaller, and can be fit with a Planck function as the time bin becomes small 
enough. This is in stark contrast to GRB 080916C, which shows no evidence of "narrowing" 
with the reducing time bin. This suggests that the Band-function component seen in GRB 
080916C is physically different from the MeV component seen in GRB 090902B. 

• We tentatively propose that phenomenologically there can be three elemental spectral com- 
ponents (Fig l23j) . namely, (I): a Band-function component (Band) that extends to a wide 
spectral regime without "narrowing" with reduced time bins, which is likely of non-thermal 
origin; (II): a quasi-thermal component (BB) that "narrows" with reducing time bins and 
that can be reduced to a blackbody (or multi-color blackbody) function; and (III): a power- 
law component (PL) that has a positive slope in uF u space and extends to very high energy 
beyond the LAT energy band. 

• Component I (Band) is the most common spectral component, which appears in 15 of 17 
GRBs. Except GRB 090926A (which may have Component III at late times), all these GRBs 
have a Band-only spectrum in the time-resolved spectral analysis. 

• Component II (BB) shows up in the time-resolved spectral analysis of GRB 090902B and 
possibly also in GRB 090510. The MeV component of these two GRBs can be fit with a power 
law with exponential cutoff (CPL). Since data demand the superposition with an additional 
PL component (Component III), the uncertainty in the spectral index of the PL component 
makes it possible to have a range of low energy photon indices for the CPL component. In 
particular, the MeV component of GRB 090902B can be adjusted to be consistent with a 
blackbody (Plank function). This is not possible for GRB 090510, whose low energy photon 
index is softer. In any case, the MeV component of GRB 090510 may be consistent with a 
multi-color blackbody. 

• Component III (PL) shows up in both GRB 090902B and the short GRB 090510, and probably 
in the late epochs of GRB 090926A as well. It has a positive slope in vF u , which suggests 
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that most energy in this component is released near or above the high energy end of the LAT 
energy band. 

• With the above three elemental emission components, one may imagine 7 possible spectral 
combinations. Most (~ 80%) of GRBs in our sample have the Band-only spectra. GRB 
090902B has the BB+PL spectra in the time resolved spectral analyses, and GRB 090510 has 
a CPL + PL spectra. Both can be considered as the superposition between Components II 
and III. GRB 090926A may have the superposition between I and III at late epochs. Other 
combinations are not identified yet with the current analysis, but some combinations (e.g. 
I+II, I+II+III) may in principle exist. 

• LAT-band emission has a delayed onset with respect to GBM-band emission in some (but 
not all) GRBs and it usually lasts much longer. In most cases (all except GRBs 090902B 
and 090510), however, the LAT and GBM photons are consistent with belonging to the 
same spectral component, suggesting a possible common origin. For bright bursts, the LAT- 
band activities usually roughly track the GBM-band activities. In the long-term, the LAT 
and GBM lightcurves sometimes (not always) show different decaying behaviors. The LAT 
lightcurves continuously decay as a power-law up to hundreds of seconds. 

• A statistical study of the spectral parameters in our sample generally confirms the previously 
found correlations between E p and luminosity, both globally in the whole sample and indi- 
vidually within each burst. We also discover preliminary rough correlations between a and f5 
(negative correlation) and between flux and a (positive correlation). Both correlations need 
confirmation from a larger sample. 

From these results, we can draw the following physical implications regarding the nature of 
GRBs. 

The Band-only spectra are inconsistent with the simplest fireball photosphere-internal-shock 
model. This is because if the Band component is non-thermal emission from the internal shock, the 
expected photosphere emission should be very bright. A natural solution is to invoke a Poynting- 
flux-dominated flow. An alternative possibility is to interpret the Band component as the photo- 
sphere emission itself. However, the following results seem to disfavor such a possibility. (1) In 
some cases (e.g. GRB 080916C), the Band-only spectrum extends to energies as high as 10s of 
GeV; (2) The low-energy photon indices in the time-resolved spectra are typically —1, much softer 
than that expected in the photosphere models; (3) There is no evidence that the Band component 
is the temporal superposition of thermal-like emission components in the Band-only sample. We 
therefore suggest that GRB 080916C and probably all Band-only GRBs may correspond to those 
GRBs whose jet composition is dominated by a Poynting flux rather than a baryonic flux (Zhang 
& Pe'er 2009; Zhang k Yan 2011). 

The existence of a bright photosphere component in GRB 090902B (see also Ryde et al. 2010; 
Pe'er et al. 2010) suggests that the composition of this GRB is likely a hot fireball without strong 
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magnetization. It is rare, but its existence nonetheless suggests that GRB outflow composition may 
be diverse. Its associated PL component is hard to interpret, but it may be from the contributions 
of multiple non-thermal spectral components (Pe'er et al. 2010). The case of GRB 090510 may 
be similar to GRB 090902B. The low-energy spectral index of the MeV component is too shallow 
to be consistent with a blackbody, but the high-latitude emission from an instantaneously ejected 
fireball (which is relevant to short GRBs) would result in a multi-color blackbody due to the angular 
superposition effect (Pe'er & Ryde 2010). 

The delayed onset of GeV emission may be simply due to one of the following two reasons: (1) 
The particle acceleration condition may be different throughout the burst. Initially, the electron 
spectral index may be steep initially (so that GeV emission is too faint to be detected), but later 
it turns to a shallower value so that GeV emission emerges above the detector sensitivity; (2) 
Initially the ejecta may be more opaque so that there was a pair-production spectral cutoff below 
the LAT band. This cutoff energy later moves to higher energies to allow LAT photons to be 
detected. Within this picture, the electron spectral index is similar throughout the burst. There 
are other models discussed in the literature to attribute GeV emission to a different origin from 
the MeV component. This is reasonable for GRB 090510 and GRB 090902B, but for most other 
GRBs this model is contrived since the GeV emission appears as the natural extension of the MeV 
Band-function to high energies. 

The GeV emission during the prompt phase is very likely not of external forward shock origin. 
This is due to the following facts: (1) In most GRBs the entire Fermi-band emission is well fit 
by a single Band component. The GeV emission is consistent with being the extension of MeV 
to high energies. (2) During the prompt phase and except for the delayed onset in some GRBs, 
the LAT-band activities in bright GRBs generally track GBM-band activities. The latter property 
is relevant even for GRB 090902B which shows clearly two components in the spectra. (3) The 
peak of GeV lightcurve coincides the second peak of GBM lightcurve for GRB 080916C. A more 
reasonable possibility is that the GeV emission during the prompt phase has an "internal" origin 
similar to its MeV counterpart. 

The origin of the long lasting GeV afterglow after the prompt emission phase (end of the 
GBM-band emission) is unclear. If it is from the external forward shock, one needs to introduce 
abnormal shock parameters, and to argue for coincidence to connect with the internal-origin early 
GeV emission to form a simple PL decay lightcurve. Alternatively, the long lasting GeV emission 
can be also of the internal origin. Future joint Fermi/ Swift observations of the early GeV/X-ray 
afterglows of some bright GRBs will help to differentiate between these possibilities. 

The two tentative correlations (a — f3 and a-flux) proposed in this paper need to be confirmed 
with a larger data sample, and their physical implications will be discussed then. 
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Table 1. The GRBs co-detected by Fermi LAT and GBM since Fermi science operation and 

until May, 2010 
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x 10~ 5 
x 10~ 5 
.02) ± 10~ 4 
x 10- 4 
x 10~ 5 
x 10~ 5 
x 10~ 5 
x 10~ 5 
x 10~ 6 
x 10~ 5 



4 84+ - 59 
4 - O4 -0.57 
1 KK+0.16 

1 ' oo -0.11 

(1,61 ±3 

8 74+ 1 - 21 

4 4 8 +0-69 
^• 4O -0.56 

> 1.07±g" 
7 14+061 

2.06±J-*| 

,+0!44 
-0.38 
(2.10 ± 

1 93+ 08 
1 - MO -0.07 
,+0.34 
-0.27 
7+0.64 
-0.51 
-1.42 
-1.26 
1+1.07 
-0.57 
R 1 c;+ 2 - 85 
D - i0 -1.81 
,+0.12 
-0.10 



7.8V 



3.68Z 
3.17j 
7.34j 
1.21 4 



1.20"* 



Spectral Type 

BAND 

BAND 

BAND 

BAND 

BAND 

BAND 

BAND 

CPL+PL 

BAND 

BB+PL 

BAND 

BAND 

BAND 

BAND 

BAND 

BAND 

BAND 



Onset Delay E u 
Y 
Y 
Y 



600 MeV 
13.2 GeV 
3 GeV 



~ 1 GeV 
~ 1 GeV 

> 100 MeV 
~ 31 GeV 
~ 30 GeV 
33.4+J;£ GeV 
~20 GeV 

> 150 MeV 
1.2 GeV 

~ 2.2 GeV 
~ 300 MeV 
~ 800 MeV 
~ 2.6 GeV 



Not e. — R eferences: (1) GRB080825C : z, Tgp - Ivan der Horst fc Connaughtonl [120081 1; (2) GRB090916 C: z -iGreirier et alj 
fcOQSh . T 90 - iGoldstein fc van der Horstl l|2008l l. E T - lAbdo et al.l ll2009al'l: (3) GR B081024B: T 90 - lAbdo et alj feoich: 
(4) G RB081215A: Tqn - IPreecel ll2008ll; (5) GRB09032 8: Ton - iGolenetskii et al] feOQSh: [6] GRB0905 10: z - iMcBreen et al.l 
bold); (7) GRB090902 B: z - Icucchiara et all ^OOgal); (8) GRB090926 A: z - fMalesani et all j2009h : (9) GRB091003: z - 
ICucchiara et alj J2009bh : (10) GRB100414A: * - Icucchiara fc Foxl d201(Jl 
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Table 2. Time-resolved and time-integrated spectral fitting parameters of 17 Fermi/LAT GRBs. 



080825C Model : Band Function 



Soq Time a ft E a K x 2 dof 

B kcV photons Q100fce y 

~™ ™ krVrmfs,,,, 



1 


0.00-6.75 


u -°'-0.04 


-2.29 ± 0.04 


135+i U 


ul -0 .007 


147.1 


154 


2 


6.75-18.1 


-0.75 ± 0.06 


„ „_+0.09 
"0-0.07 




051+° 005 


132.7 


154 


3 


18.1-25.0 




-2 17+ 017 
2,1 '-0.08 


13ll|l 


027+ ' 009 
u - uz '-0.006 


120.1 


154 


Total 


0.00-25.0 


-0.73 ± 0.03 


, .,0 + 0.04 
z,,:5 -0.03 


148 ± 9 


058 + ° 003 


265.6 


154 



080916C Model : Band Function 



Soq Time a f! E K x dof 

s kcV k*vl™'* * @ 100keV 

"1 0.00-3.70 -0.69+oqJ -2.49+',;^ 342+|^ 047+JJ ^ 9975 124 

2 3.70-9.10 -1.14 + 0.03 -2. 32+°'°° 1680+343 0.027 + 0.001 153.0 124 

3 9.10-17.0 -1.15+. ,; ,! —2.29+g°^ 9751^^ 0.016 + 0.001 125.9 124 

4 17.0-25.0 -0.99 + 0.04 -2.27+°;^ 447 leo 0.024 + 0.001 114.3 124 

5 25.0-41.0 -1.08 + 0.03 -2.49t_o ° 666+g" 0.017 + 0.001 124.2 124 

6 41.0-66.0 -1.09 + 0.04 -2.36+QQ5 696 li28 0.010 + 0.001 162.8 124 
Total 0.00-66.0 —1.05 + 0.02 -2.30 + 0.02 664+5r 0.018 + 0.001 427.5 124 



081024B Model : Band Function 



Soq 


Time 
s 


a 


P 


-^0 
keV 


K 

photons Q100fce y 
keVrrn.-* s 


X 2 dof 


1 


-0.300-0.800 


1 15 + 11 - 14 


-2.20(fixed) 


1478+jJSi" 


0.007 ± 0.001 


353.9 208 



081215A Model : Band Function 



Soq Time a fi E K x 2 dof 

s keV photons @lmke y 
r-rm r * 



1 


0.00-1.50 


-0.65 ± 0.05 


„ „_+U.14 

Z ' 2 <-0.11 

-2 16+ 010 
^■ IO -0.08 


753+^j 1 


0.059 ± 0.002 


80.0 


71 


2 


1.50-2.28 


„ r „+0.08 
U " J/ -0.07 


280l« 


n 99 ,+0.020 
v.4zc_ 017 


63.6 


61 


3 


2.28-4.93 


-0.60 ± 0.06 


-2 34+ 09 
^■• M _0.08 


178+ 2 ° 


156+ 013 

U.10D_ Q12 


66.1 


77 


4 


4.93-5.59 


-0 49+ ' 09 
w -0-08 


2 29+ 15 




0.266_ 026 


45.0 


54 


5 


5.59-8.00 


-0 72+° 16 


9 iq+013 


102+. 22 


093 + " 1,29 


47.5 


82 


Total 


0.00-8.00 


-0.71 ± 0.03 


-2 16+ 004 
z,ID -0.03 


289± 22 


110+° 005 
u - llu — 0.004 


179.9 


86 
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Table 2 — Continued 



090217 Model : Band Function 



Soq Timo a E K x 2 dof 
s kcV k^yj^L @1 00 keV 

1 0.00-7.50 -0.59 ± 0.04 -2.56+^ ^!, 365+. 0-027 ± 0.001 165.1 156~ 

2 7.50-13.1 -0.83 ± 0.05 -2.661J]"'" 47 0±5°j 0.021 ± 0.001 135.5 156 

3 13.1-19.7 -0.96 ± 0.09 -2.38t°- 22 257+ 7 5J 0.015 ± 0.002 131.1 156 

4 19.7-30.0 -0.52+°;*| -2.22+°;^ 118 j:§| - 008 to'oo3 175 ' 4 156 
Total 0.00-30.0 -0.81 ± 0.03 -2.54+°°° 418+^„ 0.015 ± 0.001 371.6 156 



090323 Model : Band Function 



Soq Time a E a K \ 2 dof 

B keV photons @100k v 

kr. Vr.m^ s 

1 5.00-14.0 -2.58t l o'i3 792+ j 3 g 0.016 ± 0.001 98.4 125 

2 14.0-25.0 -1.11 ± 0.04 — 2.54+° ^ 826±J°i 0.017 ± 0.001 127.2 125 

3 35.0-50.0 -1.08 ±0.03 -2.64+°;^ S57+® 4 , 0.018 ± 0.001 151.5 125 

4 50.0-60.0 -0.88 ±0.04 -2. 81 jig 449±| 2 0.026 ± 0.001 115.2 125 

5 60.0-135. -1.3lj:°;°j -2.62+°;^ 0.010 ± 0.001 496.7 125 

6 135.-145. -1.30 ±0.06 -2.34+°;^ 294+1* °- 017 to.001 208 ' 3 125 
Total 0.00-150. -1.22 ±0.01 -2.68+° °° 880+° n 0.012 ± 0.001 857.3 125 



090328 Model : Band Function 



Soq Time a E a K x 2 dof 

a keV photons @wokeV 

1 3.00-8.00 -0.92+JJ-'^ — 2.38+q "^q 662+.™ 0.024 ± 0.001 ~ 188.0 217 

2 12.0-20.0 -0.96 ± 0.02 -2.381°°° 727±*° 0.024 ± 0.001 199.3 217 

3 20.0-30.0 -1.15 ± 0.03 -2.30t°°° 616±g9 0.020 ± 0.001 250.7 217 
Total 0.00-30.0 -1.05 ±0.01 -2.44+°°,| 791 tso 0.018 ± 0.001 472.5 217 



090510 Model : Cut-off Powcr-Law+Power Law 



Soq Time r C PL -En ifcPL r PL K PL X dof 

S kcV photons fmkeV photons Qlfce y 



1 0.450-0.600 -0.76 ± 0.08 2688+^" 1.85+5;J5 83.7 230 

2 0.600-0.800 -0.60t°j3 4286t^|2 °- 47 ±o'.26 _1 - 73 io.07 23 - 2 tl2:3 154 ' 9 251 

3 0.800-0.900 -0.75t°° 7 777± 4 °°° °- 97 lo'.93 _1 ' 60 -o'.O7 14 - 3 t";6 52 -° 178 

4 0.900-1.00 -1.62 ±0.06 11. 5+5-3 3 8.0 134 

Total 0.450-1.00 -0.76+°°* 3624+ 7 ™ 1.06+°^* -1.66t°'°| 11.9±g 2 215.0 272 



090626 Model : Band Function 



Soq Time a E K x 2 dof 

B keV photons mi0QkeV 

r-r^n r-m ki=V r.m.*- s 

1 0.00-9.00 -0.99+°'^ -2.47t l Q 03 193+. 44 0.079 ± 0.003 340.3 186~ 

2 15.0-20.0 -1.42 ± 0.03 -2.47l° J| 391±°°> 0.040 ± 0.002 155.6 186 

3 20.0-27.0 -1.28+°'°^ — 2.58+° g| 504+°| 0.034 ± 0.001 136.5 186 

4 30.0-40.0 -1.30 ± 0.03 -2.49+° l° 444+°;! 0.025 ± 0.001 211.7 186 

— 0.06 —50 

Total 0.00-60.0 -1.40 ±0.01 -2. 62+° °* 482+H 0.025 ± 0.001 743.3 186 



090902B Model : Black Body+Power Law 



Soq 


Time 


kT (keV) 


K BB 


r PL 




x 2 




dof 






keV 


Z.39 






photons mk ^ 








1 


0.00-1.50 


75.60±; ; f 9 


38.84+JJJ^ 


-1 


.88 ± 0.02 


43.01^ 


330 


6 


264 


2 


1.50-2.25 


98.74t 3 ;| 7 


57.13t 2 ; 2 » 


-1 


„.+0.03 
' 84 -0.04 


1 -1 + 3 
,311 -4.3 


226 


3 


237 


3 


2.25-2.81 


121.20+|°° 


84.54ti; 7 7 ° 


-1 


fi1 +0.03 
■ 81 -0.04 


27 5+ 4 - 6 


217 


.5 


238 


4 


2.81-3.23 


82.52+ 3 ;f 7 


B8.00^0| 


-1 


„ n + 0.03 
' 80 -0.04 


33.6+°- 3 


199 





217 


5 


3.23-3.83 


100. 90+|™ 


69.22+ 2 ^ 


-1 


oo + 0. 03 
' 83 -0.04 


34.7i°;° 


190 


.7 


240 


6 


3.83-4.46 


86.81± 2 ;° 2 


en m +2. '20 
60 01 -2.14 


-1 


o, + 0.03 
■ oo -0.0i 


33.4+=; 7 


218 


.3 


236 


7 


4.46-4.99 


90.79t 4 ;™ 


47.82t";°| 


-1 


qo+0.03 
S3 -0.04 


38.6+S;* 


207 


4 


225 


8 


4.99-5.45 


109.50+4jj 


88. 50+IH 


-1 


„ o + 0.04 
' 82 -0.05 


3i-B±S-.S 


185 


.3 


228 


9 


5.45-5.86 


116.20+^ 


85.70+tf 3 


-1 


' 82 -0.05 




180 


5 


227 


10 


5.86-6.28 


i32.60+ 4 ; 3 2 ° 


141.20+.|;? 7 


-1 


s1 +0.04 
' 81 -0.05 




186 


.5 


233 


11 


6.28-6.61 


157 40+ 6 - 74 


155.601 7 ; 77 


-1 


o,+0.04 
ol -0.06 


38-ot*;° 


186 


.2 


228 


12 


6.61-7.19 


171 10+ 501 




-1 


„„ + 0.02 
' 8b -0.03 


87.2+ 8 7 6 3 


229 


.0 


248 


13 


7.19-7.65 


174.20+5' 3 g 


207.90t 7 ; 37 


-1 


„ 7 +0.02 
87 -0.03 


i24.stiS;i 


231 


.3 


244 


14 


7.65-8.00 


217.80+ 7 7 ;f 9 


307.ooti 2 ;™ 


-1 


.87 ± 0.02 


203.5+^.° 


223 


.0 


243 
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Table 2 — Continued 



15 


8.00-8.50 


204.80+5 5? 

— 5.48 


288.60+??? 

— 9.01 


-1 


.91 ± 0.01 


344 


"— 15.7 


319 


9 


248 


16 


8.50-9.00 


206.60+5'?! 

— 5.83 


281.00+?.?5 
— 9.16 


— 1 


.93+?' .! 
— 0.02 


375 


7 +21.5 
-19.3 


260 


2 


249 


17 


9.00-9.50 


206.20+5-?? 


270.50+5.11 

-8.91 


— 1 


.92 ± 0.01 


445 


— 18.9 


325 


6 


248 


18 


9.50-10 





135. 90+?'?? 

— 3 . 18 


209.90+5-5? 

— . 45 


— 1 


96+0.01 
— 0.02 




2 + 28.6 
— 26.0 


271 


2 


244 


19 






168.80+ 4 '55 


236.40+I'J,? 

— 7.04 


— 1 


.94 ± 0.02 


378 


4 +23.S 
-20.9 


258 


3 


244 


20 


10.5-11 


o 


195.70+6' . 3 , 

— 5 .89 


246.60+?™ 

— 8.50 


-1 


.90 ± 0.01 


352 


5 +17.7 
— 16.0 


348 


6 


247 


21 


11.0-11 




145.20+1-™ 
— 4.34 


179.10+5?? 

— 5.81 


-1 


.93 ± 0.02 


332 


2 + 20.8 
— 18.3 


278 


5 


242 


22 


11.5-12 


o 


153.10+ 4 / 43 , 


169 30+5 5? 
uo ou — 5.56 


— 1 


.92 ± 0.02 


253 


c + 18.8 
— 16.2 


241 


9 


241 


23 


12 0-12 


4 


61.07+?^ 

-2.90 


44.61+ 2 '24 


— 1 


.90 ± 0.02 


242 


Q + ] 8 A 

-15.9 


194 


7 


214 


24 


12.4-13 


2 


35.36+?-?? 

— 0.88 


31.80+°-'?!, 

— 0.90 


-1 


.92 ± 0.01 


271 


2+12.8 
' —11.9 


324 


6 


231 


25 


13.2-13 


3 


42.30 + V^ 
— 1 .59 


87.55+??, 2 

— 3.83 


-1 


.84 ± 0.03 


213 


? +27.0 
" —22.7 


141 


4 


180 


26 


13.3-13 


g 


45 32+?!? 


57 60+? ™ 

" u — 2.72 


— 1 


.87 ± 0.02 


276 


— 20.6 


175 


3 


192 


27 


13.6-13 


g 


53.27+ 2 '? 2 
— 1.94 


69.62+^-^ 


-1 


S7 +0.02 
— 0.03 


203 


7 +20.6 
' —17.3 


169 


2 


199 


28 


13.8-14 


1 


66 19 + ? 

^ — 2 . 72 


89 79 + ^'™ 


— 1 


.84 ± 0.02 


187 


o + 15.3 
°-13.8 


275 


3 


206 


29 


14. 1-14 


2 


105.70 + ^'H? 


201.80 + g gg 


— 1 


.82 ± 0.03 


169 


fi +20.2 
— 18.2 


177 


9 


204 


30 


14.2-14 


4 


120.40 + ~~ 
v -5 .70 


10. 00 


— 1 


00+O.02 

— 0. 03 


159 


q+18.7 
— 15.2 


180. 


7 


211 


31 


14.4-14 


g 


51. 74+?- 4 5, 

' —2.30 


57.16 +2 ,-SS 

" — 2. 79 


— 1 


*S6_o Q3 


186 


0+I8.S 
°— 16.2 


164 


6 


194 


32 


14.6-14 


8 


99.11 +4 -nn 
^ — 4.00 


155. 80 + 5'~ 
— 6.57 


— 1 


.85 ± 0.03 


160 


. + 19.3 
— 15.4 


173. 


6 


211 


33 


14.8-15 


Q 


71.48+ 3 ,- 30 , 
— 3.09 


115.90+t'oo 
— 5.38 


-1 


.82 ± 0.03 


149 


0+!?? 

— 15.9 


165 


7 


196 


34 


15.0-15 




102.20 + r'o2 
' — 5.26 


220.80 + ^' 7 


— 1 


.81 ± 0.03 


159 


0+21.9 
1 - 18.3 


184 


4 


202 


35 


15.1-15 


2 


102.10+ 4 4 ? 

— 4.22 


233.10 + in'i 

— 10.1 


-1 


.81 ± 0.03 


144 


6+18.9 
— 15.4 


212 


1 


199 


36 


15.2-15 




127.0 + o'lo 
— 3.73 


223.0 + I'?o 
° — 7.18 


— 1 


oc + 0.0201 
'°— 0.0234 


160 


7+! 4 -? 

— 12.5 


216 


60 


215 


37 


15.5-15 


7 


150 70+™ 
iJ "' ' u — 5.99 


254 80 + H'™ 

^ ,ou -ll .30 


— 1 


.83 ± 0.03 


120 


. + 15.9 
— 12.4 


168 


4 


221 


38 


15.7-16 


2 


59.42+*- 8 / 1 ; 
— 1.74 


63.99 + ^'!^ 

-2.12 


— 1 


.88 ± 0.02 


169 


4+; 4 -? 

— 12.4 


197 


2 


221 


39 


16.2-16 


3 


84.53 + o'?^ 

— 3.69 


132.10 + ?no 

— 6.08 


— 1 


.84 ± 0.03 


168 


q+20.3 

— 16.5 


190 


3 


203 


40 






90.82+3-67 


160.90 + «'?o 

— 6.63 


-1 


.83 ± 0.03 


158 


,+18.3 
— 15.4 


177 


3 


206 


41 


16.5-16 


7 


94.44+ 4 -55 
— 4.25 


143.00 + I'^ 

— 6.81 


-1 


.84 ± 0.03 


160 


6+19.1 
— 15.8 


169 


6 


210 


42 


16.7-16 


Q 


78.69+ 4 - 4 ^ 
' —4.10 


96.94 + ^'™ 

— ;>,29 


— 1 


oo+0. 03 
*° — 0.04 


137 


2+18.4 
— 15.1 


155 


4 


198 


43 


16.9-17 




47 97+?-?5 


40 30 + ~'~? 


— 1 


84+??? 
-0.03 


138 


7+1S-? 
— 13.1 


144 


2 


191 


44 






63 52+?-^ 


75 35+?'?? 
' ^ OJ — 2.87 


— 1 


86+??? 
— 0.03 


148 


0+15.6 
— 13.1 


171 


4 


206 


45 


17.5-17 


8 


68.97 + io« 

— ,i .26 


54.62+?'?5 

— 2 . 70 


— 1 


.85+??? 
— 0.03 


113 


7 +12.6 
— 10.6 


191 


9 


209 


46 


17.8-18. 


3 


46.21+1^ 
— 1.50 


38.75+1'?? 
°_ 1.36 


— 1 


.87 ± 0.02 


142 


O+10.4 
-9.5 


248 





228 


47 


18.3-18. 


9 




52 36+!'?? 
^ ou — 1.75 


— 1 


.88 ± 0.02 


166 


4+l° Q 6 
-9.8 


334. 





233 


48 


18.9-19. 


4 


57.29+V" 

— 1.87 


49.10+i- 7 5. 
— 1.71 


-1 


.88 ± 0.02 


156 


—9.7 


302 


1 


220 


49 


19.4-19. 


.6 


49.44+1'^ 


81.63^3-39 


-1 


.83 ± 0.03 


147 


4+18.8 


167 


7 


189 


50 


19.6-19 


7 


54.68+m 


88.95+4^ 


-1 


„,+0.02 
•°°-0.03 


164 


q+18.'9 
a -16.3 


171 


8 


192 


51 


19.7-19 


9 


57.57_ 243 


94.S9+. 4 ; 2 ? 


-1 


oc+0.02 
oa -0.03 


178 


n +i8.o 

"-16.0 


202 


2 


194 


52 


19.9-20 


1 


79 81 + 4 - 16 
' 2 - 8i -3.90 


Si-Ssts'.OS 


-1 


o, + 0.02 
oo -0.03 


197 


0+20.5 

S -17.5 


170 


6 


196 


53 


20.1-20. 


3 


43.33^.37 


42.35+. 2 ;f 8 


-1 


„ o + 0.02 
S2 -0.03 


136 


„+16.6 
°-14.9 


165 


1 


189 


54 


20.3-20 


6 


50.94tJ|2 


53.85i 2 ;« 4 


-1 


.86 ± 0.02 


193 


q+17-2 
a -15.3 


221 


4 


205 


55 


20.6-20 


9 


46.04+1;" 


Sl-23+l-il 


-1 


.87 ± 0.02 


192 


, + 16. 7 
°-14.8 


192 


6 


196 


56 


20.9-21 





42.49± 2 ; 2 4 ' 


55.46t 2 -;?° 


-1 


.84 ± 0.03 


148 


C+18.6 
a -16.0 


171 


3 


183 


57 


21.0-21 


3 


oc ,1-7+2.44 


23.88+Jg 


-1 


s _+0.02 
01 -0.03 


152 


Q+17.0 
y -14.5 


143 


5 


189 


58 


21.3-21 


7 


42 84+ 1 ' 23 


50.72+^7 


-1 


00+0.02 
• 55 -0.03 


155 


,+ 14. 8 
^-12. 7 


186 


5 


212 


59 


21.7-21 


9 


47 05+ 2 - 89 


46.i9+ 2 ;f 


-1 


84+0-02 
"*-0.03 


161 


g+17.8 
a -15.5 


152. 


6 


195 


60 


21.9-22. 


2 


49.53+^9 


42.03+. 2 ;? 4 


-1 


04+O.O2 
° -0.03 


153 


g+17.5 
D -15.1 


147 


1 


188 


61 


22.2-23 





,+4.08 


°'^-0.60 


-1 


.90 ± 0.02 


126 


0+10.2 


187 


3 


233 


Total 


0.00-30 





96 71+ - 461 
»D. <±_ 484 


71 65+ - 34 


-1 


.93 ± 0.01 


175 


,+1.2 
1 -1.3 


14732.0 


276 







090902B Model : 


Band Function - 


- Power Law 












Time 
s 


a p 


Bo 
keV 


K 

photons &1 Q 0ke y 
krVrm 2 * 


IpL 


photons @lke y 
krVrm ,±£ 


x 1 


dof 


Total 


0.00-23.0 


-0.83 ± 0.01 -3.68+0- 


12 724+ 13 
20 — 12 


0.099 ± 0.001 


, OB+1H5 
-1-»5_1. 85 


43.4 ± 1.5 


2024.3 


275 



090926A Model : Band Function 



Scq Time a f3 E K x 2 dof 

a keV photons Q100fce y 

k i= V r m 2 s 

1 0.00-2.81 -0.53+? ' 1 * -2.43+''- 1 '? 235+1? 0.106 ± 0.004 189.0 210 

— 0.03 —0.05 —15 

2 2.81-3.75 -0.48 + 0.03 -2.75+°^ 255+.j 4 0.303 +2".oio 168 - 6 196 

3 3.75-5.62 -0.57 + 0.02 -2.35 + 0.02 208 + 8 0.344 + 0.009 269.1 213 

4 5.62-7.50 -0.73 + 0.02 -2.501°;^ 326 + 15 0.191 + 0.004 229.7 210 

5 7.50-9.38 -0.63 + 0.03 -2.S,l+° ^l s 183+1 °- 255 5?o'.oo8 169 - 6 209 

6 9.38-11.2 -0.75 + 0.02 -2.52+°;J° 193+| 0.327+0;oJg 228.1 213 

7 11.2-13.1 -0.80 + 0.03 -2.29+006 154+^J 0.242+0;°i2 186.1 212 

8 13.1-15.9 -0.99 + 0.05 -2.36+. ,; 22 ISljlig 0.081+.°;°°^ 164.7 213 

9 15.9-20.0 -1.26 + 0.08 -2.07+°^ 216+6| 0.025+"™ 4 170.9 214 
Total 0.00-20.0 -0.74 + 0.01 -2.34 + 0.01 226 ± 4 0.165 + 0.002 777.1 216 



091003 Model : Band Function 
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Table 2 — Continued 



Scq 


Time 

s 


oc 





E 
keV 


K 

photons Q 100fce y 


x 2 


dof 


1 

12 
3 
4 

Total 


7.00-15.0 
15.0-18.0 
18.0-20.0 
20.0-26.0 
0.00-26.0 


-1.33 ± 0.05 
-1.01 ± 0.04 
-0.85 ± 0.03 
-1 36+0-06 

-1 09+ 02 
i ' UJ -0.01 


-2 41+ U - 2U 
„ „+0.19 
'■°^ — 0.X0 
, --+0.10 

, - s +0.05 


426+ ^ 

3sr±S 

150 ' -26 
429+ 143 

«4l 2 26 


0.012 ± 0.001 
0.040 ± 0.002 
0.094 ± 0.003 
0.014 ± 0.001 
0.024 ± 0.001 


234.5 
152.4 
218.9 
189.2 
446.2 


246 
243 
242 
246 
246 






091031 Model : Banc 


[ Function 








Scq 


Time 
s 


oc 




kcV 


K 

Photons @100keV 

kf.V (Wl^B 


x 2 


dof 


1 
2 
3 

Total 


0.00-8.00 
8.00-15.0 
15.0-25.0 
0.00-25.0 


-0.89 ± 0.06 

n oc + 0.06 
-0.Sfc._o. 05 

n 7R+ - 11 
-0.78_o 10 

„ O7+0 04 

~ U ' 8 '-0.03 


-2 44+ lu)y 

-2 50 +0 - 13 
, --+0.26 
z -°°-0.12 
, -(-+0.06 


357±|| 
467+157 

4 58l 5 33 


0.013 ± 0.001 
0.020 ± 0.001 
0.006 ± 0.001 
0.012 ± 0.001 


177.1 
173.3 
187.1 
347.2 


186 
186 
186 
186 


100116A Model : Band Function 


Seq 


Time 
s 


oc 




keV 


K 

photons @10(jkeV 
kr.V wn-i ft 


x 2 


dof 


1 

2 
3 
4 

Total 


-2.00-5.00 

80.0-90.0 

90.0-95.0 

95.0-110. 

0.00-110. 


-1 03+ Ula 

-1 03+ 05 
1 ' ua -0.04 

-1.00 ± 0.01 

-1.03 ± 0.05 
n +0.01 
^ — 0.02 


o 00+ ■ 

-2 63+ 23 
, ,,+0.11 
a ' 1J -0.09 


384+ 2 " 1 
aM -124 

791+1 92 

,al -142 

1459+m 

677llf 
2867151° 


0.006 ± 0.001 
0.010 ± 0.001 
0.033 ± 0.001 
0.009 ± 0.001 
0.004 ± 0.001 


104.8 
127.8 
156.9 
127.0 
415.6 


155 
155 
155 
155 
155 


100225A Model : Band Function 


Seq 


Time 
8 


OC 




E 
keV 


K 

photons Ql00fce y 


x 2 


dof 


1 

2 

Total 


0.00-6.00 
6.00-12.0 
0.00-12.0 


-0 53+ 11 - 22 

-0 93+ 15 
ua,:> -0.13 

-0 77+ ' 12 
u - ' '-0.11 


-2 43+U-'* 7 
^■* J -0.19 

-2 30 +0 - 26 
, , 7 +0.18 
"'-O.IO 


263+™ 

507 +3B1 
ou '-181 

375+ i29 


0.010 ± 0.002 

009+ 002 
u.uua_ 001 

0.010 ± 0.001 


51.8 
40.3 
64.5 


94 
93 
94 


100325A Model : Band Function 


Seq 


Time 
s 


oc 





Bo 
keV 


K 

photons @100keV 


x 2 


dof 


1 


-3.00-10.0 


-0 72+"- 11 


— 2 60+ i-Ky 


155 + 32 


0.014 ± 0.002 


151.6 


125 


100414A Model : Band Function 


Seq 


Time 
s 


oc 


P 


Bo 
keV 


K 

photons @100keV 


x 2 


dof 



1 1.00-7.25 -0.19lo 05 - 2 - 54 -o'w 256 -20 0.036 ± 0.002 124.3 156 

2 7.25-14.3 -0.25+°°l -2.891";^ 2 8ll2o °' 040 to'ooi 12A - b 156 

3 14.3-19.6 -0.56lo'°3 - 2 ' 53 lo'lo 361 1_6 0.047 + 0.002 135.1 156 

4 19.6-25.5 -0.76 + 0.03 -2.1$+°^ 3861™ 0.052 + 0.002 131.9 156 
Total 1.00-26.0 -0.52 + 0.02 -2.62lo"o5 3441 J 2 0.042 + 0.001 281.7 156 



Note. — Spectral models used in this paper, i.e, Band Function (BAND), Black-Body (BB), Cut-off Power Law (CPL), and Power-Law (PL), 
correspond to grbm, bbody, cutoffpl and powcrlaw in XSPEC package, respectively. Details of the formulae of these models can be found at 
http : / /heasarc . gsf c . nasa. gov/docs/xanadu/xspec/manual/XspecModels . html 
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Table 3. 


Temporal decay slopes and the time integrated photon indices of the long-term LAT 






count rate lightcurves 


Name 


f^LAT 


Tlat 


080825C 


— 0.47 ± 0.74 


-1.71 


080916C 


— 1.33 ± 0.08 


-1.77 


081024B 


-1.37 ±0.41 


-1.98 


081215A 






090217 


-0.81 ± 0.23 


-1.97 


090323 


-0.52 ± 0.67 


-1.75 


090328 


-0.96 ± 0.44 


-1.82 


090510 


-1.70 ±0.08 


-1.94 


090626 




-1.53 


090902B 


-1.40 ±0.06 


-1.76 


090926A 


-2.05 ± 0.14 


-2.03 


091003 


< -0.93 


-1.74 


091031 


-0.57 ±0.28 


-1.73 


100116A 




-1.68 


100225A 




-1.77 


100325A 


< -1.04 


-1.53 


100414A 


-1.64 ±0.89 


-1.85 
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Fig. 1. — Joint temporal and spectral analysis of GBM and LAT data for GRB 080825C. Left 
panels: the background-subtracted GBM and LAT lightcurves (from top: 8-150 keV, 150-300 keV, 
300 keV - 1 MeV, 1-30 MeV, MOO MeV), and evolution of spectra parameters (a, P, E p ). Right 
panels: an example (the brightest episode) of the observed photon spectrum as compared with the 
spectral model (top), the best fit uF u spectra of all time bins (middle), and the comparison between 
the GBM (green) and LAT (blue) count rate lightcurves in log-scale (bottom). 



-39- 



0809 16C 




20 40 60 80 Time Since Trigger (sec) 

Time Since Trigger (sec) 



Fig. 2.— Same as Figure 1, but for GRB 080916C. 
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Fig. 3.— Same as Figure 1, but for GRB 081024B. 
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Fig. 4. — Same as Figure 1, but for GRB 081215A. This burst was at an angle of 86 degrees to 
the LAT boresight. The data cannot be obtained with the standard analysis procedures. Using a 
non-standard data selection, over 100 counts above background were detected within a 0.5 s interval 
in coincidence with the main GBM peak (McEnery et al. 2008). We thus add this GRB in our 
sample, but do not add its LAT data in our analysis. 
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Fig. 5.— Same as Figure 1, but for GRB 090217. 



-43- 



-3 r 



20 8L 



090323 




Mm Uk 



13 







> 100M;V(LAT) 



an 



-4 



4- 



50 100 150 200 

Time Since Trigger (sec) 




10000.00 - 
1000.00 - 
100.00 - 
10.00 - 
1.00- 
0.10 - 
0.01 L 



10 2 10" 

Energy (keV) 



10 100 1000 

Time Since Trigger (sec) 



Fig. 6.— Same as Figure 1, but for GRB 090323. 
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Fig. 7.— Same as Figure 1, but for GRB 090328. 
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Fig. 8. — Same as Figure 1, but for GRB 090510. The applied model is cut-off power-law plus 
power-law (CPL + PL). 
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Fig. 9.— Same as Figure 1, but for GRB 090626. 
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Fig. 10. — Same as Fig. 1, but for GRB 090902B. The applied model is blackbody plus power law 
(BB + PL). 
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Fig. 11.— Same as Figure 1, but for GRB 090926A. 
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Fig. 12.— Same as Figure 1, but for GRB 091003. 
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Fig. 13.— Same as Figure 1, but for GRB 091031. 
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Fig. 14.— Same as Figure 1, but for GRB 100116A. 
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Fig. 15.— Same as Figure 1, but for GRB 100225A. 
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Fig. 16.— Same as Figure 1, but for GRB 100325A. 
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Fig. 17.— Same as Figure 1, but for GRB 100414A. 
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Fig. 18.— A comparison between GRB 080916C and GRB 090902B. Upper panel: The case of GRB 
080916C. The Band parameters are (a,/3) = (-1.0, -2. 27)(-1.0, -2.29), (-1.12, -2.25) for 0-20 s, 
2-10 s, and 3.5-8 s, respectively. Little spectral parameter variation is seen with reducing time 
bins. Lower panel: The case of GRB 090902B. (1) For 0-20 s, the Band+PL model (a = -0.58, 
P = -3.32, r PL = -2.0 with x 2 /dof = 3.52) and the CPL+PL model (r CPL = -0.59, r PL = -2.0 
with x 2 dof = 3.7) give marginally acceptable fits to the data. The CPL+PL model with Tcpl = 1 
(Rayleigh-Jeans) and the BB+PL model give unacceptable fits. (2) For 8.5-11.5 s, the Band+PL 
model (a = -0.07, (3 = -3.69, r PL = -1.97 with x 2 /dof = 1.26) and the CPL+PL model 
(r CPL = -0.08, r PL = -2.1 with x 2 dof = 1.3) give acceptable fits to the data. The CPL+PL model 
with Tc P l = 1 (x 2 /dof = 3.7) and the BB+PL model (x 2 /dof = 4.9) give marginally acceptable 
fits. (3) 9.5-10 s, the Band+PL model (a = 0.07, p < -5, r PL = -2.05 with x 2 /dof = 0.69) 
can only give an upper limit on p. The CPL+PL model (Fc P l = —0.0004, T P l = —2.1 with 
X 2 dof = 0.63) give marginally acceptable fit to the data. On the other hand, the CPL+PL model 
with r CP L = 1 (x 2 /dof = 0.92) and the BB+PL model (x 2 /dof = 1.11) give acceptable fits. Clear 
narrowing trend is seen when the time bins get smaller. 
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Fig. 19. — Distributions of the Band-function parameters a, ft, and E p in our sample (red) in 
comparison with the BATSE bright sources sample (green). The BATSE sample is adopted from 
Preece et al. (2000). 
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Fig. 20. — A comparison between the lightcurves of the blackbody component (red) and the power- 
law component (green) in GRB 090902B. The total lightcurve (the sum of the two components, 
dotted line) is also shown for comparison. 
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Fig. 21- 

correlation (panel b) for the 8 Fermi/LAT GRBs with known redshifts. The grey dots in (a) are 
previous bursts taken from Zhang et al. (2009). 





Fig. 22. — The two dimension plots of various pairs of spectral parameters, (a) a — /3, with linear 
Pearson correlation coeffcients for individual bursts marked in the inset; (b) a— flux, with linear 
Pearson correlation coeffcients for individual bursts marked in the inset; (c) E p — a; (d) E p — j3. 
For those burst without redshift, z = 2.0 is assumed (grey symbols and lines). 
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Fig. 23. — A cartoon picture of three elemental spectral components that shape GRB prompt 
emission spectra: (I) a Band-function component that is likely of the non-thermal origin; (II) a 
quasi-thermal component; and (III) an extra power-law component that extends to high energy, 
which is expected to have a cut-off near or above the high energy end of the LAT energy band. 
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Fig. 24. — Five possible spectral combinations with the three spectral components. 
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Fig. 25. — LAT photon arrival time distribution for GRB 080916C. A rough trend of gradual 
increase of the maximum photon energy with time is seen. 



